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1. 0 INTRODUCTION 


The original research contract on "Analysis of Thermal 
Stresses and Metal Movement during Welding” was initiated 
on May 15, 1969 with a appropriation of $35,390. The 
study was completed on October 14, 1970, and the final 
report of the study was published on December 15, 1970, 
as the NASA Contractor Report CR-61351. 

- fc — =^i — i 

The extension of the contract became effective on 
October 14, 1970 with an appropriation of $19,735. In 
June 1972, the contract was further extended with an 
additional appropriation of $10,000. In June 1973, the 
contract was again extended until June 30, 1974 with no 
additional funding. This report covers the work from 
October 14, 1970 to date, 
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2.0 SCOPE OF THE EXTENDED RESEARCH CONTRACT 


The scope of the extended research contract is described 
in Modification No. 3 of Contract No. NAS8-24255, issued 
in November 1970. The revised contract covered the 
following phases and tasks : 

Phase B , Task 4 . To develop a system of mathematical 
solutions and computer programs for one (1) dimensional 
analysis and for two (2) dimensional analysis. 

Phase D. Adapt Phase B system of mathematical solution 
of temperature change, stress and strain to refrac- 
tory metal alloys Cb752 and Ta222 as have been developed 
for aluminum. The emphasis of analyses should be 
made for GTWA process, 20-60 inches per minute, 
0.10-0.020 inches thick and heat input widely ranging 
about 200,000 joules/inch/inch. Necessary metal 
movement and stress measurements or materials will 
be furnished by the government as determined by MSFC. 

As decided in Modification No. 7, issued in June 1972, 
the scope of research was expanded to cover the following 
phases : 

Phase E . Perform a mathematical study of bending 
stresses a\? distortion of 2219-T87 in the thickness 
direction-aluminum only. 

Phase F. Perform a mathematical study of thermal 


stresses and metal movement in joining thin cylindri- 
cal shells-theoretical portion only. 

In June 1973, the program has been further expanded 
to include: 

Phase G . Perform experiments on thermal stresses and 
metal movement in joining thin cylindrical shells. 
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3.0 PROGRESS OF RESEARCH 

Phase B, Task 4 and Phase D 

Appendix A of this report, which is the thesis 
entitled “Investigation of Thermal Stress and Buckling 
During Welding of Tantalum and Columbium Sheet" prepared 
by Mrs. K. Anne S. Hirscn describes the work on Phase B, 

Task 4, and Phase D. 

Chapter II of the thesis describes mathematical analyses 
on thermal stresses. 

Chapters III and IV describe experimental procedures 
and results. Experimental results were compared with 
analytical predictions. Important results obtained are 
as follows: 

(1) Temperature changes during welding of tantalum 
and columbium sheets were comparable to those 
predicted by the mathematical analysis. 

(2) The one-dimensional program was found to be 
very accurate in predicting thermal strains 
parallel to the weld. Experimental results 

and theoretical results were in excellent agreement. 

(3) Since the one-dimensional program was capable 

of analyzing only one strain (or stress) component, 
a study was made to develop a two-dimensional 
program capable of analyzing all the three 
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components under the plane-stress condition. 

To overcome mathematical difficulties involved 
in the two-dimensional analysis, the finite- 
element analysis was used. The problem, however, 
is a compromise between the accuracy of calculation 
and computation cost. The accuracy of calculation 
can be improved by using a finer grid system, but 
as the number of elements increases, computation 
cost increases tremendously. The two-dimensional 
program as developed by January, 1972, was not 
able to provide results with satisfactory accuracy. 


Phase E 

Appendix B of this report, which is the thesis 
entitled "Analysis of Two-Dimensional Thermal Strains and 
Metal Movement During Welding" prepared by LCDR Jon J. Bryan, 
USN, describes the work on Phase E. 

Chapter II of the thesis describes mathematical 
formulations of the two-dimensional finite-element analysis 
of thermal stresses during welding. The significant contri- 
bution that LCDR Bryan has made is to develop a mathe- 
matical system for expressing temperature dependencies 
of various material properties, including yield stress, 
coefficient of linear thermal expansion, thermal conductivi- 
ty, etc. 
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Chapters III, IV, and V describe experimental procedures, 
data reproduction, and results, respectively. Experiments 
were conducted on butt welds and bead-on-plate welds in 
aluminum. Measurements were made of temperature changes, 
thermal strains and relative movements between both sides 
of the weld line. Computer programs were developed to 
process a large number of data generated. An important 
scope of the study by LCDR Bryan was to generate reliable 
experimental data. Important results obtained are as 
follows : 

(1) Butt vs. Bead-on-plate Characteristics . Figures 
41 through 44 of Bryan's thesis provide tempera- 
ture distribution data for the tests with strain 
gages installed and Figures 45 through 51 provide 
data on transverse and longitudinal strain dis- 
tribution resulting from these tests. Ideally, 
a butt weld becomes more and more like a bead- on- 
plate weld for multipass situations. One would 
hope for a two-dimensional model which could cover 
both cases. However, Figures 43 and 44 show a 
different initial temperature distribution for 
bead-on-plate and butt welds. This implies that 
for the plate thickness utilized in this study, 
which was 1/4 inch, the two types of weld will 
behave somewhat differently in the way of thermally- 
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induced strains, moments, and distortions. The 
bead-on-plate weld is not a two-dimensional 
situation, but the butt weld clearly is. 

This is because the arc force penetrates into 
the joint gap and the plate edges in the joint 
are evenly heated in the thickness direction. 
Figures 45, 46, 48, and 49 seem to agree with 
this observation for strains, since bead-on-plate 
values do not match the top ones . 

(2) Use of Strain Invariant . When one dea)s with 

transient thermal stresses during welding plates, 

he faces a difficult task of analysing three 

stress components, cr.., cr and t , which change 

x y xy 

with time. During the experiments at M.I.T., 
difficulties were experienced in analyzing data 
and obtaining some trends. 

In order to simplify the analysis, an attempt 
was made to utilize strain invariant, I, as 
follows (see Equation (64) in page 81) : 

T _ , 2 „ „ , 2 , . 2 . 1/2 

I = (a —oo + o + 3 t ) 

v x x y y xy ' 

On the basis of the theory of plasticity, the 

material is in the plastic condition when I is 

larger than the value of yield stress. 

Figures 56, 57, and 58 show changes of l/°y S 
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during welding. Values of I /a give rather 

ys 

consistent, explainable results compared to other 

values such as a a. and the direction of pron- 

x y 

cipal stress (6 ), shown in Figures 45 through 
55. Figures 56, 57, and 58 also show some regions 
near the weld undergo plastic deformation when 
the temperatures reach their maximum values. 

The study by LCDR Bryan provided useful experimental data. 
However, he did not have enough time to compare experimental 
data with analytical predictions. This effort is being 
carried out by Dr. T. Muraki. 

Phase F. 

Efforts have been made by Dr. T. Muraki to develop 
mathematical analysis of thermal stresses and metal move- 
ment in joining thin cylindrical shells. Basic equations 
have been developed and numerical computations are being 
made. 


T 
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4.0 PLAN FOR THE FUTURE STUDY 


A plan for the future study from July 1, 1973 through 
June 30, 1974 is as follows: 

(1) Phase G experiment 

(2) Complete description of theoretical studies 

4.1 Phase G 

Although basic equations have been developed to analyze 
thermal stresses and metal movement in joining thin cylin- 
drical shells, it is important to compare analytical 
solutions with experimental results. Efforts will be 
made to generate experimental data on cylindrical shells. 

The extent of the efforts, however, depends upon the 
availability of naval officer students. 

4.2 Complete Description of Theoretical Studies 

As described earlier in this report, a major problem 
in the theoretical study is how to develop a computer 
program with satisfactory accuracy and reasonable cost. 
During the entire course of this research, persistent 
efforts have been made to achieve that goal and computer 
progrmas have been modified for many times. 

At the end of the current research program, a complete 
description will be given of most progressed computer 
programs. 
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BY 

KATHRYN ANNE STREET HIRSCH 
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ABSTRACT 

The heat flow and stress-strain analyses of bead-on-plate 
welds is discussed and current efforts to apply the theory are 
outlined. Temperature and strain data from welding experi- 
ments on tantalum and columbium (niobium) sheet are presented. 

The experimental data are compared to analytical 
predictions obtained from both one- and two-dimensional com- 

i 

puter programs developed for the National Aeronautics and 
Space Administration. Results indicate that the one- 
dimensional program can be used to predict and analyze bead- 
on-plate weldments. The two-dimensional program, however, 
needs more programmer work in order to become a viable tool. 
Analytical and experimental evidence indicate that tantalum 
and columbium may be welded in thin sheet form but with a 
measure of local .thermal buckling. 
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I INTRODUCTION 


A. Background 

Welding is one of the most widely used methods for the 
fabrication of complex structures, because welding offers 
various advantages over other fabrication processes such as 
riveting and casting. There are, however, several undesired 
effects inherent in this process. Among these are residual 
stresses, distortion, and buckling. The primary cause of 
these phenomena is the presence of thermally induced plastic 
strain brought about during the welding process . It is the 
object of this thesis to study the effects of elevated tempera- 
ture on the amount of residual stress, distortion and buckling 
in the material. 

Residual stresses in a material may lead to low stress 
brittle fracture. Below a certain transition temperature, 
normally ductile materials may fracture catastrophically in 
brittle mode in the presence of high tensile stresses and 
notch defects. Such defects are common in welding and are 
present as cracks or lack-of-fusion spots. Since residual 
stresses in the heat affected sone may be at or near the yield 
strength of the material, very little applied stress may be 
necessary to cause fracture. 

Distortion and buckling reduce joint strength by causing 
mismatching. They may also impart initial deflection in 
structural members which can result in premature yielding. 
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Work has already been done at M.I.T. developing one-and 
two-dimensional computer analyses of strains induced by non- 
uniform heating. These programs have been compared with ex- 
perimental data obtained from aluminum and steel sheet. This 
thesis was undertaken with the thought in mind of broadening 
the previous scope by considering new materials and adding a 
buckling analysis to the wealth of strain information already 
obtained. 

Tantalum and columbium e3£hibit excel lent high temperature 
properties. Figure 1 compares the yield strength versus 
temperature curves of tantalum and columbium to those several 
of other structural materials . This is particularly important 
for aerospace applications where reentry temperatures range 
from 2000° to 4000° F depending on the trajectory. The melt- 
ing point of tantalum is 5425° F — well within the safe range. 
It exhibits a low ductile-to-brittle transition temperature 
and excellent fabricability for a refractory metal, Columbium 
(also called niobium) melts at 4474° F and is currently being 
used in several experimental high- temperature structural 
alloys. In comparison to the previously studied alloys of 
aluminum and steel, these should give investigators an idea 
of the promise of refractory metals as structural materials . 

As represented in Figure 1, the ultrahigh strength steel 
has an excellent yield stress capability but only over a 
limited temperature range. Heat treated aluminum and low 
carbon steel retain their respective strengths only as far 
along the temperature scale as does high strength steel. For 



YIELD STRESS, 10° psi 
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Figure 1 Yield Stress vs, Temperature of Several 
Structural Materials 
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any application over 1500° F these materials will suffer 
severe decreases in desirable material properties. Columbium 
and, more specifically, tantalum fall into a different cate- 
gory of materials. They demonstrate a reduction in strength 
at approximately 1500° F, however, they retain at least a 
moderate yield stress beyond this point. This indicates that 
these two materials show promise for use in applications 
which call for service in a wide range of temperature, while 
retaining structural stability. Thus, application to the 
NASA space shuttle program seems most appropriate, 

B. Welding Stress-Strain Development 

Because a weldment is locally heated by the welding arc, 
the temperature distribution in the weldment is not uniform 
and changes as welding progresses. This nonuniform tempera- 
ture distribution causes thermal stresses in the weldment 
which also change during the process. 

Figure 2 shows schematically how welding thermal 
stresses are formed. Figure 2a indicates a bead-on-plate 
weld in which a weld bead is being deposited at a speed, v. 
O-xy is the coordinate system? the origin, O, is on the sur- 
face underneath the welding arc, and the x-direction lies in 
the direction of arc travel. 

Figure 2b shows the temperature distribution along 
several cross sections. Along Section A-A, which is ahead of 
the arc, the temperature change due to welding, AT, is almost 
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zero (Figure 2b~l) . Along Section B-B, which crosses the 
welding arc, the temperature distribution is very steep 
(Figure 2b— 2 ) . Along Section C-C, some distance behind the 
arc, the distribution of temperature is as indicated in 
Figure 2b-3. Far behind the arc (Section D-D) , the tempera- 
ture change due to welding again diminishes (Figure 2b-4) . 

Figure 2c shows the distribution along these sections 

of the x-direction stress, cr j( ,. Stress in the y-direction, a , 

and shear stress, x , also exist in a two-dimensional stress 

xy 

field but are secondary and are neglected in this discussion. 
Along Section A-A, thermal stresses due to welding are essen- 
tially zero (Figure 2c-l) . The stress distribution along 
Section B-B is shown in Figure 2c-2. In the area beneath the 
welding arc stresses are near zero because molten metal cannot 
support loads. Immediately outside the weld fuddle, stresses 
are compressive because thermal expansion of these areas is 
restrained by surrounding areas that are heated to lower 
temperatures. Since the temperatures of the areas immediately 
adjacent to the puddle are quite high and the strength of the 
material is correspondingly low, stresses in these areas are 
as high as the yield strength and plastic straining occurs. 
Stresses in areas away from the weld are tensile and balance 
with the compressive stresses near the weld. In other words, 

/a ■ dy = 0, across Section B-B. 

The distribution along Section C-C is shown in Figure 2 j- 3. 
Here the weld metal and base metal regions near the weld are 
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cooling and tend to shrink. This causes tensile stresses in 
these regions near the weld and compressive stresses in areas 
further away. 

Figure 2c-4 shows the stress distribution along 
Section D-D. Continued cooling and shrinkage has left very 
high tensile stresses in and near the weld, and offsetting 
compressive stresses across the rest of the section. This 
is the residual stress distribution after complete cool-down. 

Note that the cross-hatched area, MM 1 , in Figure 2a 
indicates the region where plastic deformation occurs during 
the welding thermal cycle. The region outside MM' remains 
elastic during the entire cycle. 

C • Previous Investigations 

Under the sponsorship of NASA an ongoing series of 
projects on welding parameters has been undertaken. These 
projects comprise the first concise look at welding technology 
in this country as seen by the engineer rather than the welder 

The first project in the series consisted of an 
improvement of earlier attempts by Battelle Memorial Institute 
Columbus Laboratories to develop a computer program for cal- 
culating longitudinal stresses during bead-on-plate welds. 

The results of the Battelle study are covered in a report 
(RSIC-820) published by the Redstone Scientific Information 
Center, U. S. Army Missile Command. ^ The M.I.T. work added 
to the .. cope of the previous program by considering, in 
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addition, an investigation of tho temperature changes caused 
by the welding arc, including an analysis of temperature dis- 
tribution due to the heat generated by the welding arc. A 
system of mathematical statements was developed describing 
the phenomenon of thermal stresses and plastic strains during 
welding. Computer programs for a one-dimensional analysis of 

the problem and subsequently a finite element two-dimensional 

r 21 

analysis were written. 

The next phase of the program was to compare the 

analytic results with experimental data. Bead-on-plate welds 

were made on 2219 aluminum alloy plate, 1/4-inch thick. 

Various welding conditions were used during the experiments. 

The correlation of data proved the computer program to be an 

accurate and highly useful tool for predicting the magnitudes 

r 31 

and directions of residual stresses. In general, longitu- 

dinal strains (along the axis of the weld) were predominant. 
Transverse and shear strains were of smaller magnitude, except 
in the immediate area of the welding arc. For this reason, 
the one-dimensional computer predictions were essentially 
verified. Heat input significantly affected the extent of the 
tensile residual stress zone. High compressive stresses 
occurred in areas ahead of the moving arc. Figure 3 is a 
schematic of the development of the current welding analysis . 

The next step in the development of a suitable weld 
analysis technique was to experiment on different materials 
under different conditions. Transient strain and temperature 
data were obtained from welding experiments on low-carbon and 
high-strength steels with different strength levels (up to 
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180/000 psi yield strength). These experiments were designed 
to approximate ship structural weldments including thi c k- 
section, multi-pass butt welds. The experimental data were 
compared to the computer programs and results indicated that 
the programs could be used to analyze complex structural weld- 
ments applicable to ship and submarine fabrication. 

It was now apparent that the programs had merit for weld 
analysis, therefore, another attempt at fundamental applica- 
tion of the computer was proposed, Johnson ^ performed experi- 
ments on flame heating to remove the residual stresses induced 
by weldments in steel. He again used thick plate sections. 

The program was not as successful as previously, due to the 
large area of the flame heat source. 

D . Objective of This Study 

This particular investigation is an attempt to broaden 
the scope of previous investigations even further. Two new 
materials are being used--namely , tantalum and columbium. 

These refractory metals have previously defied all attempts 
at welding due to the ease at which they oxidixe in the pres- 
ence of elevated temperatures . Thin sheets are being used on 
the order of .012 to .015 inches thick instead of the 1/4-inch 
plate sections used in the past. Thi 3 means that the necessity 
has arisen for which we must take heat loss by convection from 
the sheet surface into consideration. All equations and the 
subsequent computer program must be modified to include these 
convection terms. 
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This ifivuy ligation is of particular inioroaL I < j HAiiA 
because they are considering manufacturing the space shuttle 
from tantalum or columbium. In sponsoring the contract, they 
requested the experimental welds be made using a gas tungsten 
arc apparatus in order that this report might simulate welds 
which could be made on a production-line basis. 

As a follow up to previous studies using the NASA 
computer program, this investigation will use the one- and 
two-dimensional weld analyses. 

Some study of buckling of sheet after welding has 
already been done ^ but little progress was made. This 
paper devotes a section to an attempt to develop an extremely 
simplified buckling analysis applicable to welded structures. 

As far as we know, this is the first study ever published 
covering transient local buckling near the arc. Studies on 
transient thermal strains near the arc exist as do works on 

r 6 1 

buckling of the entire plate after welding, L J however, none 
has been published on local buckling which moves with the arc. 
This is a significant problem in the fabrication of structures 
using thin sheet. Without the use of a computer, such a study 
is virtually impossible. 
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II THEORY 


A. Heat rlow Analysis 


Derivations of the heat flow equations have been set 

[ 7 ] 

forth by several authors. Chief among these are Myers, et al 
who give a derivation and a number of solutions which are 
applicable to various welding processes. Unfortunately, it is 
not yet possible to calculate with satisfactory accuracy, the 
thermal cycles in the vicinity of the weld made under given 
welding conditions . 

Figure 4 shows schematically the temperature distribution 
in a plate on whose surface a weld bead is being laid at a 
speed, v. Curves 1 to 6 represent isothermal curves on the 
surface, while the dotted curves represent isothermal curves 
on the transverse section, ABCD. O-xy is the coordinate axis; 
the origin, 0, is on the surface underneath the welding arc, 
the x-axis lies in the direction of welding, and the z-axis is 
placed in the thickness of the plate, downward. 

The fundamental expression for heat conduction is given 
by the Fourier heat flow equation, as follows: 


30 

3t 


o 

3^0 


3x‘ 


2 2 
3^0 , 3^0 

* — o — 2 

dy* 


The mathematical analysis of heat flow in a weldment is 
essentially a solution of this equation for a given initial 
condition {initial temperature distribution) and a boundary 
condition {shape and intensity of the heat source, geometry 
of the weldment, etc.). 
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Tw o features of heat flow during metal-arc welding are: 

1) The heat source moves, usually at a constant speed, 
on or near the surface of the workpiece. 

2) The size of the heat source (welding arc) is small 
compared to the size of the workpiece. 


Consider the case of a solid through which heat is flowing 
but in which no heat is generated. The temperature 0 at the 
point P(x,y,z) will be a continuous function of x, y, z and t, 
and the first differential coefficients of G will also be 
continuous . If we define Q = dQ/dt to be the rate of heat 
transfer across one face of a rectangular parallelepiped, 
where p is the density and c the specific heat of the solid, 
we have the time rate of change of internal energy per unit 
volume equal to: 

6 V = p ■ c ( ||) (i) 

If we expand this expression: 


,99. „ 9 ■ 3 /\9G\ , 9 

C + 3? (: V 


( 2 ) 


This then becomes 

„ ^ ,90\ - w 3 2 0 . 9 2 0 . 9 2 0x . 9X 
P * C (-rr) = X ( 5* + g- + g-) + Ko- 

3x 2 3y 2 3z 2 36 


[« 


> ! * #>" * <H' 1 ’ 


(3) 


We assume that 3X/99 = 0, thus removing the last terms of the 
equation. 

Such a condition is realized in many cases of practical 
importance . 


The temperature distribution for welding is calculated 
from the analytic solution of the linearized heat flow equation 
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’or a moving line source. This equation is linear because 
the material properties have been assumed to be constant with 
temperature. The linear heat flow equation for a flat plate 
is : 


ae 

at 


cs 




a 2 e 


dz‘ 


( 4 ) 


where 


0 = temperature change 
k = X/c , thermal diffusivity 

P 

c = specific heat 
p = density 

X = thermal conductivity 


Quasi-stationary . In the case of a moving heat source, 
it is convenient to express this equation in a coordinate 
system moving with the heat source . The equation then 


becomes : 


2 2 2 

3 Z 9 ^ 3 4 G . a 0 

— 2 + "~2 + ' — 2 

3w^ dy 3z 4 


v 3_9_ 

tc 3w 


( 5 ) 


where 


w = x - vt 

V = speed or moving source 

If this equation is solved for a moving line source of 
intensity, q, where 
q « Q/h 

H = plate thickness 
Q = net heat input 
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The following solution is obtained 

_v 

6 - ®o “ jh e 


2K 


w 


K o 


( 6 ) 


in which 

v - /n? 2 + y 2 

k q {z) is the modified Bessel function of the second kind 
and zero order. 


Nonstationary . For two-dimensional nonstationary heat 
flow, consider an instantaneous line source, q? cal/cm occur- 
ring at P' (x 1 ,y') in an infinite plate at a time, t' , and 
then extinguishing the equation for temperature change at 
P(x,y) and at time, T, then becomes! 

_ (x-x 1 ) 2 + (y-y 1 ) 2 

di 

(7) 


6 - 8 o - 




4itk (t-t ’ ) 


Cp 


Temperature change due to a moving point source can now 
be obtained using equation (7) in the integrated form: 


t „ (x - vt) 2 + y 2 
o e 4k {t + t, - t) 

e - e = f 2 i 

° 0 4ktt (t + t, - t) 


■“3. 

cp 


( 8 ) 


Equation (8) can also be expressed as follows: 

..2 


_V *i< ( W 

e - fl . JL e " 2K / 

o 4irX 0 

Yl t 

4k r l 


- S 


_ i i 

% 


dK 


<9> 


where w - x - v (t -ft,) 

o 1 

Y 2 = (v/4k) 2 (w 2 + y 2 ) 
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Equation (9) represents a general solution for 
two-dimensional heat flow when laying a straight weld bead. 

By assuming that welding continued from t = 0 to t Q = + 60 / 
the temperature in the nonstationary state is obtained as 
follows : 


0 



v . 

2k* J 



- -*w 

3 e 2l< k (~ r) 
2ttX e k o K 2k } 


( 10 ) 


Consider a plate lying in the xy-plane and its thickness, 
b in the direction of z to be so small that the temperature 
may be taken to be constant over it. This is a plausible as- 
sumption using the generalized plane problem: 


! b/2 

4T / 0 (x,y , z) dz = 6(x,y) (11) 

D -b/2 


Let E be -the outer conductivity of tlv j material, X its thermal 
conductivity, p its density, and c its specific heat, then 
the differential equation satisfied by the temperature in the 
plate is found to be: 


3 2 0 9 2 0 

2 2 
3jT 3y^ 


pc 3 0 2E 

TT Tt “ 1“ 



0 


(12) 


where q o is the temperature of the surrounding medium. 

2 f 8 1 

Defining a , the radiation constant, as : L J 

2 2E 
pcb 


( 13 ) 
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The equation becomes: 



24 - a 2 e - 


0 


(14) 


For the sake of simplicity# we will use 9 instead of 0 in the 
rest of this thesis. 

If we solve this new radiation equation using the modified 
Bessel function as we did before, the new solution becomes: 


v 

— • *' 1 * 

0 " 6 o “ jfl e ZK <o (/ST 75? r) 


where 


2 a * , e , 2 v 

T ~ 4 ^ (w + y ) 


a . ^ a 2 + v_ 
a a * Ik 


(15) 


a 2 is defined as the radiation constant with units of 1 /sec 
and is obtained from the formula: 
a 2 = 2E/cpb 


where 


E ” radiation heat from the unit surface by unit 

2 

temperature difference per unit time, cal/ cm sec 
b = thickness of the plate, cm. 

E can be determined experimentally from the formula 

c 


[93 


w 

E - s 


0 - e 




( 16 ) 


where 

w = weight of specimen# gm 
s - surface area# cm 2 
0 = temperature of plate 
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0 q - room temperature 
V c = cooling rate, 0 c/sec 

By this method we can treat the heat flow in a thin sheet as 

2 

a two-dimensional heat flow problem with a representing a heat 

2 

sink. In this case, a is not constant but varies with 
temperature . 

The problem is not purely heat loss due to radiation into 
the atmosphere, because copper backing plates were used in the 
escperiment. These undoubtedly had an effect on heat dissipa- 
tion. Since the problem is somewhat complicated for a precise 

2 [ 17 ] 

analysis, we simply use a to represent the heat sink, 

B. Stress Analysis 

The second major aspect of the analysis of welding 
processes is the calculation of stresses and strains which 
result from temperature distributions. 

Although many problems may involve only elastic effects, 
considerable work has been done on plastic deformation. In 
general, metals behave elastically up to some limiting value 
of stress or equivalently some limiting value of strain. For 
most, this behavior may be considered to be linear, that is, 
the stress is proportional to the strain* The constant of 
proportionality is the Young's modulus of the material. If 
the material is strained beyond the yield point, the stress 
increases, but at a lower rate than before. For some materials 
very little strain occurs, and they may be considered perfectly 
plastic. 
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If the stress is then removed, the material behaves 
elastically and returns to the Kero stress level with a finite 
amount of strain. Wh^n the load is again applied, the mate- 
rial behaves elastically up to the maximum previous stress. 

When a material has undergone plastic deformation, this 
deformation remains after the load has been removed. Although 
the phenomenon is not time dependent, the material does remem- 
ber its past history. It is possible to obtain different dis- 
tributions of strain for the same stress distributions if 
loading histories are different. As a result, stress and 
strain must be calculated for small increments of loading not 
merely for the final loading if plastic deformation occurs. 

There are several different theories for predicting the 
yield point of a material which has first undergone plastic 
deformation in tension and is then loaded in compression or 
vice versa. They are illustrated in Figure 5. For a purely 
plastic material, all of these theories are equivalent. The 
present analysis employs the second theory in which the 
initial yield points are independent. 

A normal stress-strain curv..: is obtained from a uniaxial 
tension test. If such tests are run at different temperatures, 
different curves will result. In general, both the yield 
strength and Young's modulus will decrease with temperature. 

In this way a family of curves is obtained which may be 
pictured as a three-dimensional surface with the third coor- 
dinate representing temperature as shown in Figure 6. This 
surface, however, has been constructed by considering the 
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material behavior at a fixed temperature only, and does not 
really indicate how the material behaves under changing 
temperatures , 

As has been mentioned in the section on heat flow, the 
two-dimensional approximation gives a fairly realistic model 
of the welding of thin sheet since the temperature does not 
vary with thickness. In this case the plane stress formula- 
tion may be used with stress and strain being considered 

uniform in the thickness direction, and a , the stress in 

z 

this direction assumed to be zero. 

The equilibrium equations involve the longitudinal and 

transverse stresses rt and a and the shear stress t , these 

x y xy 

equations are : 


3a 


3x 

3a 


x j 
— + 


3t 

Si ~ C 

3y 


(17) 

(18) 


3t 

+ -JSL 

3y 3x 

where there are no body forces acting on the material . 

In order for the material to remain continuous, the 
longitudinal and transverse strains and e , and the shear 
stra.' 


e , must satisfy the relation: 


3 2 c. 


x 


2 2 
3 e 3 e 

X _ 2 2K. =, o 


3y- 


3x* 


3x3y 


compatibility equation (19) 


The stresses and strains are related by the following: 


= 1/E(a x - ya ) + dT + ej 


+ AeP 
x 


( 20 ) 
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e y *= 1/E (o - ya x ) + aT + + Ae^ 


'xy 


1/G t -f + Ae^ 

' 1 scy 3ty xy 


( 21 ) 

( 22 ) 


where 


E = young 1 s modulus 
y = Poisson’s ratio 
G - shear modulus - E/2 (1 + y) 
a = coefficient of linear thermal expansion 
T ~ temperature above reference temperature 

Ae^# Ae^, Ae^ are the plastic strains due to the last 
x y 3cy 

increment of load » 

e* 3 , e^, e? are the plastic strains due to previous 

x y 3cy 

increments of load. 


The incremented plastic strains are found from the 
stresses using the following relations s 
P Ae 

ie x = it Ua x " °y > 
e J 

Ae 

Ae^ — n? — ^ {2c ~ c ) 

y 2c y u x' 


( 23 ) 


( 24 ) 


Ae P ■ 3 T 

xy 2 o T xy 


( 25 ) 


where 


a e " 


t 2 2 , o 2 , 

(a x + a y " a x a y + 3T xy ) 


1/2 


( 26 ) 


, — 2 9 9 1/2 

4£ p - 2//r + Ca e 5> + ieE Ac? + (AcL) ] 


X 


( 27 ) 
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There are several types of strains involved in this 
discussion. The first type is the total strains/ c , etc. 

These strains are related to the actual displacements in the 
material and are the strains that are actually measured. 

The thermal strain, aT, is related to the expansion that 
would take place if each part of the material were completely 
free to expand. That is, the surfaces of the material are not 
restrained, and the temperature distribution is such that the 
thermal strains satisfy the compatibility equations. 

For most temperature distributions, however, the thermal 
strains alone do not satisfy compatibility, and additional 
strain is present. This is the mechanical strain or stress- 
producing strain. This strain has, in general, both an elas- 
tic and a plastic component and is related to stress by the 
stress-strain diagram and the Prandtl-Reuss relation. These 
stresses must satisfy the equilibrium equations and boundary 
conditions involving surface forces. 

The problem of an infinitely long plate with width, 2c, 
has been solved in the following manner. Using the stress- 
strain relation: 

£ = 1/E a + aT + (28) 

X X X 

if we introduce the nondimensional quantities: 

S - a /o e=e/e t = aT/E e„ = e£/ 

x o x o ' o p x o 

n = y/c H = E/E Q (29) 

where a is the yield stress at some reference temperature, 
o 

e Q r and, E q is the Young's modulus at that temperature, the 
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stress-strain relation becomes: 


e x - S/H + t + e^ 


( 30 ) 


The integrals for the net force and the moments become 
1 

/ S * dy » 0 (31) 

-1 

1 

f S ■ y • dy « 0 (32) 

-1 


If we consider 
bead-on-plate weld, 
follows : 
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Therefore, we can assume 

°x = f( y> 

and from the equilibrium condition 


3ct 
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tjX 
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(36) 


(37) 


Then 

s2e x 

— » o (38) 

3y 

This is considered the one-dimensional analysis of welding. 
The theory is not applicable in the region near the welding 
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arc, however, it can be used as an approximation for a stress 
field in which 

3cf v 3<y v 

”5x << ~5y < 39 > 

i 


This analysis was developed by Tall and was the forerunner 
of the original research program at Battel le Memorial 
Institute. 

As previously stated, the compatibility equation ^ 
reduces to 


which implies that the strain distribution is linear. That 
is: e x * a + by (40) 

When the linear form of the strain is substituted into 
the integral equations using the stress-strain relation, the 
resulting equations may be solved for the unknown coefficients 
a and b. 

In this way, the final expression for the total strain 
is obtained as: 

1 

e x (A 1 “ A 2 n)I { H(T + 

1 

, “ < A 2 “ A 3 n) 1 £ H(T + e p>^ dTl1 

where 

1 2 

/ Hn ■ dn 

A , m — 

X 1 1,1, 

/ Hdn / Hn dn - ( / Hndn) 

-i -l -l 


(41) 


{ 42 ) 
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11 * 1 
/ Hdn / Hn dn - ( / Hgdn) 
-1 -1 -1 


l 

/ Hdn 



/ Hdn f Hn dn - ( / Hndn) 

-l -i -1 


(44) 


The nondimensional thermal strain, x, and the variation 
of the nondimensional Young's modulus, H, are calculated from 
the temperature distribution, T. Since the plastic strain, 

E , is not known, the strain cannot be calculated directly; 
however, an initial approximation to the total strain can be 
obtained by assuming that the plastic strain is zero. The 
mechanical strain is then calculated from the expression 

e n = e x - T (45) 

This mechanical strain may be used to obtain a first 
approximation to the plastic strain using the stress-strain 
curve as shown in Figure 7 . 

C , Stress Calculation 

At this point it would perhaps, be helpful to understand 
exactly how the computer uses the stress-strain calculations 
to produce results. 

The time units and corresponding temperature are 
inputted into the program and may come from the analytic 



38 


Figure 7 
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Obtaining Plastic Strain from Stress-Strain Curve 



temperature program or experimental data. The time at which 
the arc reaches the gages is also read in. The material 
properties used include Young’s modulas, E, the initial yield 
stress, cr f the coefficient of linear thermal expansion, a, 
and a strain-hardening parameter, m. 

The program performs calculations at several 
nondimensional transverse positions. These positions range 
from the weld center line to the edge of the plate with 
finer divisions nearer the center line. 

In the expression for total strain, e , the integrals 
are taken across the entire plate, however, the program uses 
values across only half the plate. The resulting expression 
for total strain ahead of the weld is given by: 


= (A 1 -A 3 n) / H(T+e )dTi + (A 4 *.i-A 2 ) / H (T+e ) dy (46) 


where 


/ Hr) (n - t) dri 
0 
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(47) 
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/ H (n - 1/2) dr] 
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f Hdn 
0 
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(49) 


(50) 


where 


/ Hdrt I Hn (n-1/2) dri - / Hndn / H(n-l/2)dn 
0 0 0 0 


D 


(51) 
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Behind the weld the total strain is given by the simpler 
expression: 

1 

C = A / H(T + c )dn (52) 

* Q P 

where 

A = r A— (53) 

/ Hdn 
0 

In the program these integrals are evaluated using Simpson’s 
Rule . 

At the start of the calculation, the material is assumed 
to be free of any plastic strain. The plastic strain is 
accumulated during successive iterations until reverse yield- 
ing occurs . 

At each time step the plastic strain is calculated using 
the mechanical strain component only of the total measured 
amount. Figure 8 illustrates several possible situations of 
elastic and plastic strains as well as reverse yielding. For 
clarity the curves are shown as being independent of tempera- 
ture but the actual calculation takes such changes into 
account. 

Each value of plastic strain calculated is compared with 
the previous one until either subsequent values agree within 
one percent or the program has performed twenty iterations 
and ceases. 

If convergence is obtained, the stress for each point 
is calculated and printed out along with total, mechanical, 
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and plastic strains. When the last time step is completed, 
the calculation is repeated at all points at the reference 
temperature t0° F) . This gives the residual stress 
distribution. 

Two-dimensional Program * The current two-dimer, sional 
program considers the average stress over an entire block 
in a preset grid. As the arc is moved closer to that par- 
ticular grid block, the stress is calculated. The grid 
shape has a finer spacing nearer the weld to give more pre- 
cise values in the area of greatest interest. One major 
disadvantage of a finite element program is the cost- 
effective decision necessary for determining the fineness of 
the grid. Obviously, for best results the grid should be 
made as small as possible* This is, unfortunately, not 
possible where money for running programs is limited? there- 
fore some compromise must be made. 

First, the temperature distribution around the moving 
arc is calculated. Then the stress field is divided into 
a set of similar blocks of width, h Qr shown in Figure 9. 
The time intervals represented by the block width must be 
short enough so that the temperature and thermal stress for 
each increment may be regarded as being constant. Since the 
greatest changes in temperature occur near the arc, narrow 
strips are used in areas near the arc. 

The calculation starts on a block some distance ahead of 
the welding arc where the temperature change is negligible. 
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and the stresses are purely elastic. Time zero is fixed on 
the block* 

First stresses in the block crossing the origin, O, are 
calculated based on elasticity theory. Then stresses in the 
second block are calculated by adding stresses, due to the 
temperature increment. In this case, analysis is made whether 
or not any plastic deformation takes place. It is assumed 
that the amount of stresses at a given point does not exceed 
the yield stress of the material at the temperature of that 
point. Similar analyses are conducted step by step on the 
following blocks. Thus the stress distribution in the entire 
field is determined. ^ 11 ^ 

i). Thermal Buckling 

/ 

Buckling type distortion of welded plate can be 
classified into the following two types. First is the buck- 
ling of the whole plate after welding has been completed. 

The second is local buckling near the arc during welding. 

In practice we clamp the plate to avoid buckling. There- 
fore, in most applications the first type is the problem. 

The plate does not buckle while it is clamped, but when it 
is released. This phenomenon has been studied by both 
Masuhuchi and Watanabe and SatohP'^llo work has been done on 
the second type,* because it requires complicated analysis. 

.. « it 

In the experiments conducted in this study, buckling of 
the entire plate was not observed, probably because the small 
weld produced by the gas tungsten arc process did not produce 
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enough residual compressive stresses. Local buckling was 
observed. Figure 10 illustrates its appearance. 

The complete analysis of local buckling is very 
complicated. The ideal solution would be to develop a simple 
model to determine the critical free distance, that is, the 
distance between the clamps. Only a simple introduction 
is developed here. 

In welding of a thin plate, buckling may occur due to 
compressive thermal stresses during cooling. Compressive 
thermal stress is considered to be caused by inherent shrink- 
age due to welding. As a result, the critical inherent 
shrinkage of buckling is obtained by assuming that compressive 
thermal stress is equivalent to stress in the plate produced 
by concentrated compressive forces, P - Eh(e^), acting at 
both ends of the weld line. 

Where 

E * Young's Modulus 
h «= plate thickness (cm) 

= inherent shrinkage (cm) 

When the value of P_„ (critical) is reached, buckling occurs. 

Ci 

Stresses can be obtained using the theory of elasticity. For 
an aspect ratio (length/width) less than or equal to 2, the 
equation is given by: 

a v (at the reference points) *= 

oo 

g^-tl 4-4 I (-l) n/ ^ (mru 4- l)e nlTV ] (54) 

n-2,4,6 


where V - I»/B. 




Figure 10 Buckling Appearance after Welding 
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The formula for inherent shrinkage is: 




kirV 

12(1 - v 2 )B 


(55) 


where 

v = Poisson’s ratio 

k » numerical factor dependent upon the aspect ratio L/B. 


k is approximately represented by: 
k ** 11 + 50(B/L) 2 


(56) 


and is independent of material. 

It is considered that a plate will buckle in the 
longitudinal direction when the magnitude of exceeds the 
critical value determined by the previous equation. Watanabe 
and Satoh analyzed the relationship between inherent shrinkage 
and welding conditions in carbon steel. They found the 
formula to be: 


U-l) 


* 0.136 x 10 


-6 


I 
h/v 


■v2 


(57) 


where 

X = welding current (A) 
v = arc travel speed (cn\/sec) 

Extensive experimental evidence is necessary to determine a 
similar formula for other materials. 

Consider a long strip of width 2a, length, L, and welded 
along the center line, (bee Figure 11) The following 
assumptions are made regarding residual stress components: 
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-O* |y| > c 

longitudinal stress, 0 ^ » « 

* o" 0 < | y | < c 


transverse, « 0 
shear, t , * 0 

Since the residual stress must be symmetric: 
a 

• / crdy * 0 
0 x 

and 
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(59) 

(60) 

(61) 

(62) 


If we denote the deformation in the z-direct..\on as w, the 
equation of equilibrium is: 


D = 


/9 4 w , , 

( i ? +2 


+ ii*> « * afw 

1 


— 5 
3x^3y 


where 


D « rigidity of the plate 


Eh' 


12(1 - v*) 


(63) 


(64) 


also called bending stiffness. 


°x h 


For boundary conditions, it was assumed that the plate is 
simply supported at both ends, x *= 0, L and that the plate 
is free along y « ±A. 

If we consider a plate symmetrically heated and supported 
on rigid edge members, the most important temperature para- 
meters are the average temperature rise of the plate over that 
of the edge member, and the nonuniformity parameter representing 
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the difference between the maximum and minimum temperature 
rise. 

A temperature distribution of the form 

T (y> “ T o “ iT< b ) (65> 

is chosen to represent the tent-like distribution obtained 
in welding. 

The boundary conditions for a clamped edge determine both 
deflection and slope normal to the boundary to be zero: 

a) - 0 || s || cos a + |£ sin a - 0 (66) 


where a is the angle between the normal, n, and the direction 
of x. (See Figure 12 . ) 

Thus, the bending moments m x , and m^, and the twisting 
moment m^ are given by: 

m x = my - (1 +v)DoM 0 (67) 
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where 

m 6 
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f 0zdz 
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temperature moment 


In most pure buckling problems m 0 * 0. 
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Hence, 


( 68 ) 


(69) 


also 


(70) 

(71) 


( 72 ) 
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where 


+h/2 

n e B E / edz 

-h/2 


Mean temperature 


(73) 


6 a temperature distribution 

e e °, e ® ** strain components in the midplane 
xx yy xy 


n x K DEe xx + ve yy “ (1 + v ) an e J 
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If the plate is held fixed along its entire edge 
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XX 




(77) 


n„ = n„ « -(1 + v)akn Q (78) 

x y y 

V = o <”> 

where 

k -s — — k -— stretching stiffness (80) 

1 - 

Since in thermal buckling problems in contrast to 

ordinary buckling problems, the membrane forces n„, n , and 

x y 

n are not constants, but functions of x, y — expecially in 
xy 

this case of the .tent-like temperature distribution — it is, in 
general, not possible to find exact solutions. The added bur- 
den of the clamped boundary conditions puts this problem well 
into the realm of a doctorate thesis by itself and. 
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subsequently will not be attempted here. There is, 
unfortunately, no reference available today which even attempts 
to solve this problem. Parts, however, can be found in 
Parkus, Johns, Timoshenko, Nowacki, Cox. 13,14 ' 15,16 ^ 
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III PROCEDURES 
A. Strain Measurement 

Therroal strains in welding are measured frequently by the 
use of electrical resistance strain gages . These operate on 
the principle that certain conductors exhibit a change in 
electrical resistance with a change in strain. Gages are 
mounted on the test specimens and the resistance variation 
across the gage is measured as welding progresses. In the 
case of thermal strains in welding# the observed resistance 
change# AR, consists of: 

AR « AR x (e e ) + AR 2 (s p ) + AR 3 («T) + AR^ (T ) 

where 

ARi (e^) = the resistance change corresponding to elastic 
mechanical strain# e e . 

AR 2 (e p ) = the resistance change corresponding to plastic 
mechanical strain, 

P 

AR^taT) = the resistance change corresponding to 
temperature induced thermal strain# aT. 

AR^ (T) « the resistance change caused by thermo-electric 

effects in the gage itself. 

It is not possible to separate strain into its elastic 
and plastic components, but AR 3 (aT) and Ar^ (t) can be subtrac- 
ted by either calculation or experiment. This so-called "ap- 
parent strain" may be measured by heating a mounted strain 
gage in an oven and obtaining a plot of apparent strain versus 
temperature . 
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A schematic diagram of the experimental set-up is shown 
in Figure 13. One strain gage and thermocouple was cut out 
from each of the original tantalum and niobium sheets to which 
a weld bead had been applied. The pieces cut out were appro- 
ximately two inches by three inches in size so that the strain 
recorded by each of the gages would have a negligible mechani- 
cal component. The pieces of tantalum and niobium with their 
respective gages, along with a third unmounted strain gage and 
thermocouple were then placed into a furnace. As the tempera- 
ture in the furnace was raised from room temperature to 1000 °F, 
strain change (resistance change) of each of the electric- 
resistance strain gages was measured via a Wheatstone bridge 
set up. The induced EMF created by the temperature gradient 
between the thermocouples in the furnace and a 32° F reference 
bath was also measured, from which temperature was determined. 

The data plotted as temperature versus strain is shown 
in Figure 14. 

B. Apparatus and Procedure 

1. Test specimens . Plate dimensions and gage locations 
are shown in Figure 15. The arrangement of constraining 
clamps is shown in Figure 16. The room temperature material 
properties of the test plates are found in Table 1. 

2. Strain gages and thermocouples . The gages used in 
this investigation were type HT-1212-5B, manufactured by BLH 
Electronics, Waltham, Massachusetts. These are high tempera- 
ture, free-filament gages which include, a thermocouple as an 




Figure 13 Schematic of 
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Figure 14 Temperature vs Strain curve for Thermocouple 
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Figure 15 Specimen Dimensions 
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Figure 16 clamp Arrangement 


TABLE 1 


Material 

Tensile 

0.2% 

Young 1 s 

Coefficient 

Density 

Thermal 

Specific 


Strength 

psi 

Yield 

Strength 

psi 

Modulus 

psi 

of Linear 
Thermal 
Expansion 
per °C 

g/cm 3 

Conductivity 

cal/sec/craVoc/cm 

Heat 

cal/gm °C 

Tantalum 

55,000 

- 

27 x 10 6 

6.5 x 1Q“ 6 

16.6 

0.13 

0.034 

Columbium 

85,000 

32,000 

22.7 x io^ 

7.06 x 10“ 6 

8.57 

0.13 

0.067 


cn 

o 


integral part of the mechanism. Gage properties are; 
Designation HT-1212-5B 

Gird Length 5/16 inch 


Grid Width 


3/32 inch 


Temperature Range 
Resistance 
Gage Factor 
Cement 


-320° to +1200° F 
120 ohms 
4.11 ±1% 

Rokide — BLH 


3. Instrumentation . Strain gages were connected into a 
potentiometric circuit (half-Wheatstone bridge) and calibrated 
as indicated in Figure 17. Figure 18 shows the calibration 
circuit of the thermocouples which were referenced to a 32° F 
ice bath. Both circuits were fed into a Honeywell continuous 
recording 12-channel Visicorder. Temperature and strain were 
recorded simultaneously at the gage location. Weld passes were 
timed by means of an electric stop watch as well as a timer 
integral to the chart recorder. 


C. Welding Equipment and Conditions 

The power source used was an ac-dc Heliwelder (manual 
GTA) manufactured by the Air Reduction Company , The torch 
from this machine was wired to an automatic beam and carriage 
electronic governor manufactured by the Linde Division of 
Union Carbide Corporation. This was done to insure constant 
travel speed which could not be guaranteed using only the 
manual apparatus. The travel speed, arc voltage, and amperage 


Figure 17 Strain Gage ciicuit 


STRAIN 

GAGE 


DUMMY 

RESISTANCE 



cr> 

to 





Figure 19 Thermocouple Circuit 


WHEELCO POTENTIOMETER 
MODEL# 320P SN.06JI487 



' REFERENCE 
JUNCTION — 32°F 




Cf% 

to 
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were preset before each pass. Arc length was adjustable 
during the weld. 

D, Experimental Procedures 

The experimental operation is shown schematically in 
Figure 19. The test plate, i, «■ rumented at BLH, was clamped 
between copper plates to insure no buckling during the welding 
process. Welding speed was preset at 8 ipm, arc voltage at 10, 
and amperage at 40. The visicorder was actuated and the arc 
was struck. As the welding torch began to move down the 
plate, the arc length was adjusted and the timer started. The 
recorder output was marked when the arc passed the strain gage 
location. When the welding head reached the end of the plate 
the arc was extinguished and the plate allowed to cool . The 
recorder continued to monitor the gages for approximately 
three minutes until conditions appeared stable. Periodical 
readings were taken until the sheet cooled to room temperature. 
The clamps were then released while the recorder was still 
monitoring. 
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Figure 19 Schematic of Apparatus and Procedure 
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IV EXPERIMENTAL RESULTS 

The results of this investigation consist of mechanical 
strain data recorded during the welding of test specimens. 

The data is presented below in the form of plots of mechanical 
strain versus time and temperature versus time. For compari- 
son purposes, each experimental strain curve is preceded by 
theoretical strain curves obtained from the NASA one- 
dimensional and two-dimensional computer programs. Only longi- 
tudinal strains are plotted from the two-dimensional program 
due to the fact that the transverse strains are of a signifi- 
cantly smaller magnitude. 

The horizontal axis of all graphs is a scale of time 
measured in seconds. Zero time is arbitrary, occurring some 
time after the arc has stabilized and is moving down the plate 
toward the gage location. The point at which the arc passes 
this location is indicated on the graph. The limiting time on 
the scales merely represents a time by which the plates had 
reached steady conditions. 

The vertical scale of the first two graphs is temperature 
in degrees Fahrenheit measured from the thermocouples on the 
samples. The subsequent graphs are calibrated in p in/in 
strain on the vertical axis. Positive values indicate ten- 
sile strain and negative values indicate compressive strain. 

A. Discussion and Analysis— General 

The weld appearance showed the presence of a white 
powder which was thought .to be an oxide of the base plate. 
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When the plate was forcibly bowed after the clamps were 
removed, the weld area demonstrated severe brittleness and 
cracked in many places. The weld itself was not as continuous 
as it might have been due to the difficulty in maintaining a 
constant arc length. Several times during welding the arc 
extinguished but this is thought to have only minor effect. 

Local buckling was observed in the weld region but no 
attempt was made to measure the amount. This would have been 
impossible in this experiment due to the presence of longitu- 
dinal oxide-induced cracks along the crown of the weld bead 
which would prevent any precise measurement of deformation. 

A major purpose of this thesis was to compare experimen- 
tal results on thin sheet to the previously tested one-dimen- 
sional computer analysis and the new two-dimensional model. 

Tantalum strain gage number 3 which was located 4.5 inch 
from the weld showed a continuous tension response. The pro- 
gram predicted compression. Such disagreement did not occur 
on the columbium specimen. This can be attributed to the 
fact that tantalum gage number 3 was not located far enough 
from the weld to fall into the compression region. (See 
Figure 2.) The program had been previously used on aluminum 
and steel and used the same criterion to predict the location 
of the compressive region in tantalum. This is evidently not 
a correct assumption due to the unique temperature properties 
of this material discussed in the -previous section. 

In both tantalum and columbium gage number 3 was located 
so far from the weld center line than the copper backing 
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became a major factor and significantly lower readings were 
recorded. These should generally be ignored when analyzing 
the worth of the computer model compared with experimental 
data. 

B. Heat Flow 

In order to begin a heat flow calculation, it is 
necessary to determine a quantity known as arc efficiency, ri . 
This factor is related to heat in the following manner: 

Q = n*v.l 

where 

Q » thermal power of heat source in watts 
v = welding speed in in/sec 
I « currents in amperes 

Figures 20 through 32 show calculated temperature changes 
during welding. Measured temperature changes are shown in 
Figures 33 and 34. 

In order to determine the correct value for arc 
efficiency, n, calculations were made with three values of 
n, 0.2, 0.3, 0.4. Figures 20 through 34 show the calculations 
for each of these values . For n *= 0 . 4 the theoretical values 
were significantly higher than the experimental ones. For 
r\ = 0.3 the values were closer but not quite low enough. 

The n « 0.2 value showed closest agreement. 

For most GTA welding processes the arc efficiency is 
0.2 to 0.5. (See Figure 35.) This discrepancy was due to 
not including heat loss by radiation in‘ the analysis. 
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Calculated Temperature Change during welding 
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Figure 23 Calculated Teirsperature Change during Welding of Tantalum 
(efficiency ™ 0*1, 2.1 inch from weld center) 
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Figure 24 
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Figure 26 Calculated Temperature Change during Welding 
b*f Tantalum (efficiency ~ 0.1, 4.5 inch from 
weld center) 

Note: Time and Temperature Scales Are Reversed 
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Figure 27 Calculated Temperature Change during Welding of Columbium 
(efficiency ■=* 0.1, 0.9 inch from weld center) 



Figure 28 Calculated Temper at ur e change during V7elding of Columbium 
(efficiency - 0.2, 0,9 inch from weld center) 
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Figure 30 Calculated Temperature change during Welding of columbium 
(efficiency = 0.1, 2.1 inch from weld center) 

Note: Add 75° F to all values. 
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Figure 32 Calculated Temperature change during welding of 
Columbium (efficiency - 0.1, 4,5 inch from weld 

center) 

Note: Time and Temperature Scsles are Reversed. 


TEMPERATURE (°F) 



MID j POINT TIME (SECONDS) 

Figure 33 Experimental Temperature Change During Welding of Tantalum 
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Figure 35 Measured Values of Efficiency for Various 
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Next, the radiation term, E , was added to the temperature 
calculation program. The best curve fit obtained by varying E 
and n is given in Figure 36 . This corresponded to an E value., 

r 8 1 

of approximately 40 times that obtained by Naka 1 1 (E « watts/ 

2 

in • °F) , and an n of 0.3. This was thought to be the best 
temperature correlation possible. The figures show that the 
computed and measured temperature tends to peak slightly ear" 
lier than the calculation and decreases more rapidly. This is 
probably a further effect of heat loss from the surface. 

When the arc efficiency calculation was done at a point 
further away from the weld center line than the original cal- 
culation, the efficiency tended to take on a slightly lower 
value. This was due to the presence of the copper backing 
plate acting as a heat sink. 

C. Strains — Stress Analysis 
One-Dimensional 

Figures 37 and 38 show that longitudinal strains 
calculated by the one-dimensional computer program coincide 
extremely well with measured results (Figures 39 and 40) . 

Peaks occur approximately 5 seconds sooner after the arc passes 
in the measured graphs than they do in the calculated graphs . 
Peak heights agree within a maximum of 400 |i in/in for the 
tantalum specimens and 100 y in/in for the columbium specimens. 

The tantalum strain maximum was of a lower magnitude 
than the columbium values . This was not the trend predicted 


TEMPERATURE t°F) 



/x. INCH / INCH 


2600 

2400 

2200 

2000 

1800 

1600 

1500 

1400 

1200 

!000 

800 

600 

400 

200 


-2 00 

-400 
— 600 


120 40 

15 


80 100 120 140 160 180 200 

TIME (SECONDS) 


Figure 37 Calculated Strain Change during "elding of Tantalum 
Using the One-Dimensional Program 
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Figure 40 Experimental strain Change during Welding of Columbium 
Note; strain Scale incorrect — Read correct values 

from the numbers located on the Graph Itself* 
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in Figure 41. The error occurred due to the uncertain value 
of yield stress of columbium. For annealed products the value 
is 45,000 psi but for wrought products a high of 55,000 psi is 
reached. Figure 41 utilizes the lower value while the pro- 
gram uses the higher one . 

The peak values obtained for both specimens of 
approximately 3,000 y in/in corresponds to the values achieved 
for 1/4 inch thick 2219 aluminum alloy by Arita. This demon- 
strates the conjecture made earlier that columbium and tan- 
talum will be competitive with aluminum for structural 
purposes . 

Two-Dimensional 

In the two-dimensional calculation all strains followed 
the sarnie pattern of peak shape and all showed a similar lag 
time of 5 seconds behind measured values. The two-dimensional 
calculation results appear in Figures 42 through 45. 

There were two great disparities in the comparison of 
two-dimensional theoretical and measured results . The first 
was the actual numerical value the peak obtained. This was 
significantly {1000 y in/in compared to 3000 y in/in) lower 
than the measured value . Since the one -dimensional values 
agree with the experimental results, and taking into consi- 
deration the fact that the two-dimensional program is still 
in the developmental stage, these values can be discounted as 
programmer error. This does not negate the worth of the two- 
dimensional analysis, however, because the peak shapes and 
locations were correct. 
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Figure 42 Calculated Strain Change during welding of 
Tantalum using the Two-dimensional Program 
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Figure 44 calculated Strain Change during Welding of 
Columbium Using Twp-aimensional Program 
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Figure 45 Calculated Strain Change during Welding of 

Colurotoium Using the Two-dimensional program at 
Time Step after that used in Figure 44. 
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Also, the ratio of tantalum^^ to columbium max strain 
valurs were computed. Although the numbers did not agree, 
the peak height ratios fell in the 1.25-1.7 range the same 
as the one-dimensional and experimental results. 

The second discrepancy with experimental results occurred 
as a dip in the theoretical peak at the point where the arc 
passed the strain gages. This was due to the nature of the 
finite element program. The time steps were unfortunately 
arranged such that the computer passed from one grid block to 
another just as the torch passed the gages. This is a cor- 
rectable factor if the size and placement of grid blocks is 
changed. By ignoring the dip, the correct peak shape is 
obtained. 

D. Residual Stress 

Graphs of residual stress versus distance from the weld 
center line at successive times are shown in Figures 46 
through 50. These values were taken from the one-dimensional 
program. The fact that the area near the weld becomes tensile 
as time progresses while the area away from the weld is com- 
pressive demonstrates the validity of the discussion on resi- 
dual stress in the introductory section. The narrow width of 
the tensile zone indicates that there will be little post-weld 
buckling by the Watanabe-Satoh analysis. ^ 

Figure 51 shows how the stress field changes with time 
at three different distances from the weld. The closest 
distance of .9 inches shows the greatest variation and 
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Figure 46 Calculated Change is Residual Stress with Distance from 

Weld Center during Welding of Tantalum, 5 seconds before 
arc reaches point of measurement . 
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Figure 47 Calculated change :in Residual stress with Distance from 

Weld center during Welding of Tantalum Just as Arc Passes point of Measurement 
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Figure 48 Calculated change in Residual Stress with Distance from 
Weld Center during Welding of Tantalum, 5 seconds after 
Arc Passed Point of lleasurement 
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Figure 49 calculated change in R Q sidual Stress with Distance from 

Weld Center during Welding of Tantalum, 15 Seconds after 
Arc Passed Point of Measurement 
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Figure 50 calculated Change in Residual Stress with Distance from 

Weld Center during Welding of Tantalum, 35 seconds after 
Arc Fassed Point of Measurement 
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Figure 51 Calculated change in Residual Stress with Time 
during Welding of Tantalum 
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stabilizes in the tensile region. The farthest distance from 
the weld varies only slightly from tensile to compressive. 

This is merely another method of demonstrating the same infor- 
mation that appears in Figures 46 through 50. 
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CONCLUSIONS 

The purpose of this investigation was to determine the 
reactions produced by welding Tantalum and Columbium sheet 
and compare them with computer predictions. 

1. The temperature distribution during welding of 
Tantalum and Columbium is comparable to that expected in other 
structural materials of similar thickness. This results in a 
thermal stress of suitable proportion to allow the use of 
Tantalum and Columbium as a structural material. 

2. The weldability of Tantalum and Columbium was not 
sufficiently studied to be able to recommend the "best" tech- 
nique, however, it appears that weldability is possible using 
present, but carefully controlled, production techniques. 

3. Local thermal buckling may well be a major considera- 
tion in any attempt at fabrication using large Tantalum or 
Columbium sheet. Theoretical predictional ability is not 
currently state-of-the-art. 

4. The previously developed one-dimensional computer 
program is an excellent tool for predicting expected stresses 
and strains. Experimental results and theoretical results 
were in complete agreement. 

5. The two-dimensional program appears to be able to be 
developed in the- future into a viable prediction tool. Cur- 
rently the program needs some adjustment to give correct 
strain magnitudes. 
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RECOMMENDATIONS 

1. Future welding of tantalum and columbium should be 
done in a more controlled atmosphere to prevent oxide 
contamination . 

2* Asbestos lining should be used from the standpoint 
of analysis between any necessary backing plates and the base 
plate to prevent any unnecessary heat flow complications. 

3. An extensive theoretical analysis of local buckling 
should be undertaken in an attempt to supply a viable mathe- 
matical model which can then be experimentally tested. 

4. Further research should be done on how small the 
finite element grid size can be made and still remain cost 
effective. 

5. The premise of the two-dimensional program should 
be changed in that the arc should never have to jump from 
one block to another , but the block should move along with 


the arc. 
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erature dependent properties; transient strain and temper- 
ature distribution data obtained from experiments on 
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time approximating ship structural weldments. 

The transient strain response was found to be pre- 
dominantly longitudinal but transverse strains were signifi- 
cant in the region of the welding arc. It was also found 
that residual strains correlated reasonably between similar 
experiments with a variation of the second stress deviator 
tensor invariant J 2 representing the "state of stress" 
rather than by individual strain components. 

Several recommendations are made concerning continued 
experimental investigation aimed at further development 
of the National Aeronautics and Space Adminstration programs. 
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I : INTRODUCTION 


A. Background t 

The phenomenon of thermal stresses and strains resulting 
from welding processes lead to difficulties in the fabrication 
of structures. Non-uniform temperature distribution leads 
to residual stresses which in turn provide problems from 
stress corrosion# buckling# brittle fracture, fatigue strength, 
and distortion. Depending upon the type of material being 
welded, any or all of the above problems may be important 
considerations in the structure's design and its service per- 
formance (1) . All contribute to the reliability of the 
basic structure (2, pp. 2). Masubuchi (2 - 5), Klein and 
Masubuchi (6), and Klein (7) provide interpretive reports 
which review the particulars of each of the above phenomena 
for aluminum and steel as well as comprehensive reviews of 
the literature. Almost all of the previous studies refer 
to residual stresses and distortion and only recently has 
much attention been given to the actual thermal stresses 
during welding. This is due to the complexity of the problem 
characterised by large temperature changes in small areas 
near the welding arc with its resulting non-elastic deforma- 
tion, temperature dependency of the material properties and/or 
phase transformations, and complex boundary conditions 
resulting from conditions of geometry and multi-pass welding 
(3, pp* 2,3). 
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P Thermal Stresses During Welding , Residual Stresses and 
Distortion ; 

In order to study the effects of residual stresses 
and distortion performance of welded structures# it is first 
necessary to understand the relative magnitude and distribution 
of residual stresses and distortion in weldments. The 
mechanism by which thermal stresses and strains are developed 
in welded plates is best described by Masubuchi (8# 3, pp. 4-6) 
and is repeated here: 

Figure 1 shows schematically changes of temperature 
and stresses during welding. A bead-on-plate weld is being 
made along the x-axis. The welding arc, which is moving at 
a speed, v, is presently located at the origin, 0, as shown 
in Figure la. 

Figure lb shows temperature distribution along several 
cross section. Along Section A-A, which is ahead of the 
welding arc, the temperature change due to welding, AT, is 
almost zero. Along Section B-B, which crosses the welding 
arc, the temperature distribution is very steep. Along 
Section C-C, which is some distance behind the welding arc, 
the distribution of temperature change is as shown in 
Figure lb-3. Along Section D-D, which is very far from the 

welding arc, the temperature change due to welding again 

* 

diminishes. 

Figure lc shows the distribution of stresses along 
these sections in the x-direction, o . Stress in the 
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y-direction, a , and shearing stress, x , also exists 
y *y 

in a two-dimensional stress field** 

Along Section A-A, thermal stresses due to welding are 
almost zero. The stress distribution along Section B-B is 
shown in Figure lc-2. Stresses in regions underneath the 
welding arc are close to zero, because molten metal does 
not support loads. Stresses in regions somewhat away from 
the arc are compressive, because the expansion of these 
areas is restrained by surrounding metal that is at lower 
temperatures. Since the temperatures of these areas are 
quite high and the yield strength of the material is low, 
stresses in these areas are as high as the yield strength 
of material at corresponding temperatures. The magnitude 
of compressive stress passes through a maximum with increas- 
ing distance from the weld or with decreasing temperature. 
However, stresses in areas away from the weld are tensile 
and balance with compressive stresses in areas near the 
weld. In other words, 

fa x * dY = 0 ( 1 ) 

across Section B-B.** Thus, the stress distribution along 
* — — 

In general three-dimensional stress field, six stress 
components, a x * <j z , x xy * x 2y # x zx exist. 

Equation (1) neglects the effect of <7 and j on the 
equilibrium condition. y xy 
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Section B-B is as shown in Figure lc-2. 

Stresses are distributed along Section C-C as shown in 
Figure lc-3. Since the weld metal and base metal regions 
near the weld have cooled, they try to shrink causing 
tensile stresses in regions close to the weld. As the 
distance from the weld increases, the stresses first change 
to compressive and then become tensile. 

Figure lc-4 shows the stress distribution along Section 
D-D. High tensile stresses are produced in regions near 
the weld, while compressive stresses are produced in regions 
away from the weld. The distribution of residual stresses 
that remain after welding is. completed as is shown in the 
figure. 

The cross-hatched area, M-M, in Figure la shows the 
region where plastic deformation occurs during the welding 
thermal cycle. The ellipse near the origin 0 indicates the 
region where the metal is melted. The region outside the 
cross-hatched area remains elastic during the entire weld- 
ing thermal cycle. 

As shown in Figure 1, thermal stresses during welding 
are produced by a complex mechanism which involves plastic 
deformation at a wide range of temperatures from room temp- 
erature up to the melting temperature. Because of the 
difficulty in analyzing plastic deformation, especially at 
elevated temperatures, mathematical analyses were limited 









for very simple cases such as spot welding. 

However, on the basis of recent developments in com- 
puter technology, it appears that a technological break- 
through as far as the analysis of thermal stresses during 
welding are concerned is not too far away. 

C. Previous investigations : 

The closest approach to a breakthrough to date has 
been the one-dimensional strip theory approach first pre- 
sented by Tall (9) and extended by Masubuchi and associates 
at Battelle Memorial Institute (10) and Massachusetts 
Institute of Technology (7, 11, 12) under the support of 
NASA's George C. Marshall Space Plight Center. These 
studies have been reviewed in detail in references (3- 6) 
and will not be repeated here. 

A great deal of effort in the last several years has 
been devoted to developing finite element programs which 
will extend one- dimensional analysis to two and three 
dimensions (13 - 20) . Nearly all hav' met some qualitative 
success but in general, quantitative results have been 
lacking. The finite element method utilizes an approximate 
model to derive a set of equations which is then solved 
exactly. The flexibility of the finite element method will 
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allow the extension of the flat plate problem into more 
complex and real life geometries. It provides a complete 
strain picture , including longitudinal# transverse, and 
shear strains. It has the capability to model the restraint 
conditions 5n the actual structure and weld. Its greatest 
drawback is high expense from considerable computer time 
for large degree of freedom problems and the post analysis 
of the vast quantities of data the technique generates. 
Interactive graphics will greatly aid this latter problem 
in the future. 

Masubuchi and Iwaki (13) and Masuhuchi and Andrews (14) 
coupled the elastic-plastic analysis developed by Wang (15) 
with the thermal loading calculations for moving heat 
sources previously utilized by the one-dimensional programs. 
Reference (13) was qualitative and reference (14) provided 
some quantitative results with course meshes but appeared to 
have transition element instability when a smaller mesh was 
utilized. In this situation, an element which reaches 
yield stress is considered a transition element. However, 
in the next iteration, the equivalent stress from that 
element drops below yield and the element becomes elastic 
again. The whole situation jumps between the two sets of 
values. In order to overcome this difficulty, a scheme 
was devised to quarter the time step in which this insta- 
bility occurred. This worked for a large number of cases 



but not all cases. Sometimes increasing the time step 
also eliminated the instability. Reference (10) may offer 
a way to correct this problem. Yamada (17) seems to have 
met with more success dealing with thermal stresses due to 
rapid heating by controlling the number of elements which 
yield. 

Hibbit and Marcal (18) offer the most analytic of the 
approaches to date. Their model treats the weld process 
as a thermo-mechanical problem. A finite element formula- 
tion derived from the uncoupled thermal and mechanical 
energy balances forms the basis of the model. Tt attempts 
to deal with material non-linearity due to temperature 
dependence of thermal properties and in the fusion problem 
where material phase change xs accompanied by a latent heat 
effect. It includes radiation as a cooldown mechanism and 
finite strain effects during isothermal loading. However, 
there was little agreement with experimentally measured 
residual stresses and they concluded that their finite 
element model did not include significant material behavior 

Shinn (19) approaches the two-dimensional distortion 
of a panel structure due to welding with the assumption 
of elastic deformation during the welding process. Com- 
putations were carried out using the one- dimensional 
experimental values of unconstrained angular changes along 
the welded edge and its equivalent constrained welding 
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moment as an input to the computer program. The results 
were not completely successful because of the elastic 
assumption and the questionable accuracy of the experimental 
data input into the program, it is of the opinion of this 
author that this technique, on the basis of current tech- 
nology, will provide the most fruitful results in the 
immediate future in terms of financial and manpower invest- 
ments. The basic welding problem is highly non-linear 
and may be too complex to solve in a completely analytical 
form. A purely empirical approach does not appear very 
fruitful, either. It is not a simple task to determine 
thermal stresses in small regions heated to high temperatures. 
Without proper analysis, it would be difficult to adequately 
interpret experimental data. 

The empirical-analytic technique, however, is only as 
good as experimental data provided. At present, material 
and physical properties at the range of temperature from 
room to melting, even for the most common of structural 
materials, are sorely lacking in the literature. More 
effort must be expended in providing this critical information. 
The approximations utilized in the interim will be the 
controlling feature of this computational technique. 
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D. Rim and Purpose of the Present Study : 

Muraki and Masubuchi are presently developing a new 
elasto-plastic finite element computer program at Massachusetts 
Institute of Technology along the lines described by Yamada 
(17) with temperature distribution adopted from reference (20) . 
It is the purpose of this report to provide experimental 
data on a heat treatable tempered aluminum alloy to verify 
this program. Electric resistance thermocouples, three 
element strain gage rosetts, single element strain gages, 
and an extensioneter will be utilized to obtain transient 
and residual strain and displacement data in the center and 
near the edge of 30x18 inch panels for an automatic GMA 
welding process. Experiments will be performed for bead-on- 
plate and butt welding processes under similar heat inputs 
to verify: 

1. that experiments may be duplicated under reasonably 
similar conditions to give repeatable results; 

2. to compare butt and bead-on-plate welds for similar- 
ities and differences; 

to compare strains in both top and bottom of the 
plates to study the importance of bending strains as com- 
pared with transverse and longitudinal strains; 

4* to observe transient principal strain magnitudes 
and directions ; 

5. to measure residual stress/strain distributions 



in the longitudinal and transverse directions; 

6. to measure the deflection of the two half plates 
in a butt welding process under tack welded conditions; 

7. to present longitudinal and transverse strains for 
computational technique verification, and; 

8. to collect and present in usable form available 
physical and mechanical temperature dependent data for 
use in the computer analysis. 

Standard strain gages will be utilized which may not 
exceed 400°F. Since the material properties in the heat 
affected zone (HAZ) are not fully known or understood and 
the high temperature strain gages have proven unreliable 
(6, 7) in part because of temperature' compensation data and 
in part by the lack of material properties data within the 
HAZ, strain measurements will be beyond the HAZ. 
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II: MATERIAL BEHAVIOR 


A. General : 

A new finite element program to study residual stresses 
and metal movement is presently being developed by M.I.T.’s 
Department of Ocean Engineering concurrent with this thesis. 
Since an attempt is being made here to provide material 
properties for and consistent with this new program, a 
brief and simple formulation of the tangent stiffness 
technique, incremental stress strain relationship, and 
thermal loading terms is provided for the sole purpose of 
demonstrating the importance of temperature dependency. Then, 
material properties consistent with this formulation are 
presented for 6061 aluminum alloy in the T6 and T651 condi- 
tion, the material utilized for the experiments presented 
in Chapter III. 

B. Mathematical Formulation: 

The most common finite element approach to plasticity 
problems in welding has been the Tangent Stiffness Approach. 
The problem is first solved elastically for the general 
loads. Then the elastic limit loads are determined by 
scaling down the given loads until the element with the 
maximum equivalent stress is just reached. The differences 
between the given loads and the elastic limit loads are then 
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% 

Figure 2: Equivalent Stress and Strain (14,pp 20) 



Figure 3: Iteration Process to Bring Back the Stress Point 

on to the curve with H** The element in the tran- 

sitlon region starts out at point B as elastic* The elastic 
increment. BC of the element for the prescribed load increment 
dP is first calculated. The stress point corresponding to 
C should be 0 on the stress-strain curve* Using the modulus 
corresponding to BD, the Increment BF for dP Is next calcu- 
lated* This process may be repeated until the solution 
convergee (17* pp 300,303)* 
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divided into a predetermined number of loading increments. 

The problem is then solved including plasticity for each 
element. When an element becomes plastic, the tangent 
stiffness approach is used to determine an effective elasto- 
plastic stiffness for that element. In this approach, 
the Prandtl-Reuss relation and the plastic modulus are 
combined with the elastic stress-strain relations to produce 
an incremental stress-strain relation which includes the 
effects of plasticity. These incremental relations depend 
upon the instantaneous stress levels in each element (14, 
pp. 4-5). 

In general, 

da * E. de (2) 

where E is the tangent modulus coefficient. E. has an 

."Jr fiC 

elastic term and a plastic term. The plastic term is zero 
if no yielding takes place. Then, E fc may be used to estab- 
lish the tangent stiffness array K for each element, 

K =» I B* E. B d (Volume) (3) 

j , S rS 

Volume 

where defines the strain displacement relationships. 

Thus, where loading is applied in increments, the structure 
stiffness array may be found at each stage of loading by 
evaluating each element's tangent stiffness (21, pp. 12 - 13). 

The incremental pseudo-thermal .loading term likewise 



becomes 


dP 


n 


f de, d (Volume) 

Vofume^ * ' 


(4) 


where de^ is the thermal strain caused by the welding 
process. Summed over the entire structure , one obtains: 


K dU 
dU < 
de ' 
da 


s 

K* 1 dP 


B fc dU 


E t B dU 


<5) 

( 6 ) 

(7) 

( 8 ) 

To determine the tangent modulus coefficient related 
to equation (2) , one must consider the following relations 

Prandtl-Reuss Flow Rule: 


de p - I dip 


yielding Surface Behavior! 

•• 9a ® 
dcr ® y cIcj 


Material Behavior: 


(9) 


( 10 ) 


*0 + E t r p 


"o + H’r c 


da - E? de„ 
t p 


H‘ de. 


(ID 


Utilizing equations (9) * (12) and recalling that 

The p and e superscripts refer to plastic and elastic 
respectively* 
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G 


da » Dde * D (ds - de*) (13) 

where de is the total incremental strain and 9 is given in 




Equation (24) for the plane stress situation, one arrives 
at the following relationship: 


da 


D-D 


D 

* 

3c T 

~ 8a 

/V 

. H 1 

+ 15 

3a 


D |£ 
Z, da 




) 


(14) 


This reduces to equation (2) where the quantity in £ ^ 
is the tangent modulus coefficient. This is consistent 
with reference (17, pp. 310) . 

If we assume plane stress conditions and consider the 
equivalent stress for Von Mises* isotropic material, 

S 2 - a 2 - o x <j y + o y 2 + 3r xy 2 (15) 

Differentiating, 

. r 

da * — I (2a x - a ) da x + (2a 

2a *■ 1 1 


This may be rewritten in matrix form as 

r i 


a ) da + 6t dt I 
x' y xy xyl 

(16) 


da *= — a T 


1 0 
4 ? o 

0 0 3 


da 


(17) 


or 


da = - a T P da 


(18) 
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where 


*y 


r 


da. 


x 


(19) 


da * 1 do 




dr. 


xy 


( 20 ) 


Comparing equations (10) and (18) , we see that 


35 T 


Bo 

Transposing, 

3a 

3a 


IT 
— a f 


( 21 ) 


1 „T 1 „ 

— Fa « —Fa 

a ** ** c ** *** 


( 22 ) 


Substituting the previous relationships into equation (14) 
we arrive at 

r 

.T 


da 


D F a a F D 

D — A/ «*»> 

~ a 2 H' + a T F D F a 

AO A» 


de 


(23) 


for a two dimensional, plane stress, isotropic case, with 
the Von Mises' yield criterion where 


D * 


1 - v‘ 


F = 

A* 


1 ~2 


1 

'7 


l 

o 


1 V 
0 1 

G 0 
1 

0 

0 
3 


0 

0 

1-v 

2 


(24) 


(25) 


_ 2 2 9 5 

c = a - o a + o +3t 

x x y y xy 


( 15 ) 
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Equation (23) is but a simple example of the more 
complex general case. The important point to note here is 
that E, a, and H* are functions of temperature and this 
makes quite complex even the most simple of welding problem 
examples . 

In the previous study at the Massachusetts Institute 
of Technology, Andrews (14, pp. 25) assumed the simplest 
model possible for the nodal force matrix, i.e., that of 
a constant temperature of each element, in this case, the 
nodal force matrix for each increment becomes 

dP_ * f (B T „ E. de,) d (Volume) (4) 

f'V-P J n 

Volume 



where 

AT * average temperature change of element 
a = coefficient of linear thermal expansion 

If the element in question is plastic, the tangent stiffness 
matrix is used. Masubuchi and Iwaki (13, pp. 9) assumed 
a much more complex loading function which accounted for 
the equivalent nodal forces due to the dependency of the 
coefficient of thermal expansion and Young 1 s modulus on 
temperature. The plastic terms additionally took the 
temperature dependency of the yield surface into account. 
Yielding is complex. It can occur during the heating stage 
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and also in the subsequent cooling stage since some 
elements can be stressed to their yield point in the 
reverse direction. Yamada (17, pp. 298) assumed the same 
simplified type of expression as equation (4) in his 
study of thermal stresses due to rapid heating. As 
in the case of the incremental stress-strain relation, the 
important point to note is that in varying degrees of 
difficulty, the loading term is dependent upon E, a, and 
a which are in turn dependent upon temperature. 

The temperature distribution, uncoupled from the 
mechanical problem just described, may be approximated by 
classical, finite difference, or finite element techniques 
(9-11, 13, 14, 18, 20, 24-30). In terms of the conductiv- 
ity tensor, k, the general representation of the aniso- 
tropic heterogeneous continuum is 

dT 

pc = V • (k * VT) + U'" (25) 

where p is the density, c is the specific heat, k is the 
thermal conductivity, and U'" is the heat generation per 
unit volume (23, pp. 29, 44-45). Fourier's law applied 
to the boundary conditions provides 

q = 


-K * VT 


(26) 
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where q is the heat flux normal to the boundary surface. 

If the continuum is isotropic, these vector equations 
reduce to the familiar 

p c * V • (kVT) + U'” 

q *= -kVT 

Equation (27) may be reduced to 

pc || = Vk • VT + kV 2 T + U"' (29) 

If k is a function of space only, equation (29) is 
linear. On the other hand, where k depends on temperature 
alone, 


(27) 

(28) 



equation (29) may be modified to 


pc 


dT 

dt 


dU 

dT 


(VT) 


kv 2 t 


u 


III 


(30) 


which is nonlinear. For homogeneous isotropic continua k 
is constant and equation (29) reduces to 


dT 

dt 


aV 2 T + 


UMI 

pc 


(31) 


where a is the thermal diffusivity. Since p, c, and k 
are functions of temperature, the welding problem becomes 
highly nonlinear. Coupled with phase changes, this 
problem becomes even more complex. 


C. Material Properties : 

As discussed in Section B, tabulation of the material 
properties as a function of temperature is of prime impor- 
tance. This information is difficult to obtain and often 
simply is not available. In this experimental work, alloy 
6061 in the T651 Condition is utilized. It is one of the 
most versatile of the wrought, heat treatable aluminum 
alloys. An understanding of the behavior of this common, 
relatively inexpensive -material during welding and heat 
treating may also contribute to the understanding of the 
behavior of quenched and tempered steels such as HY-80 and 
HY-130 currently in use today. 6061 is readily welded by 
all methods and has excellent weldability characteristics. 



The plate used in this experiment has a Federal Specifi- 
cation of QQA-250/11. 

1. Composition (31# Code 3206# pp. 1) 



Percent 
Min. Max. 

Chromium 

0.15 

0.35 

Copper 

0.15 

0-40 

Iron 

- 

0.7 

Magnesium 

0.8 

1.2 

Manganese 

- 

0.15 

Silicon 

0.40 

0.8 

Titanium 

- 

0.15 

Zinc 

- 

0.25 

Other impurities 



each 


0.05 

total 

- 

0.15 

Aluminum 

Balance 


TABLE 1: Composition of 6061 Aluminum Alloy 

2. Heat Treatment (31, Code 3206# pp. 1) 

Annealed (0 Condition): Heat to 775 Q F for 2 to 6 

hours# cool at 50°F per hour maximum to 500°F. The rate 
of subsequent cooling is unimportant. 

T4 Condition : Solution heat treat to 970°F# water 

quench# and naturally age at room temperature to a sub- 
stantially stable condition. 

T6 Condition : Solution heat treat to 970°F# water 

quench# artificially age by precipitation heat treatment 
at 320°F.for 16 to 20 hours or 350°F 6 to 10 hours and 
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air cool. 

T 651/T451 Condition Same as the T6/T4 Condition followed 
by a 1 1/2 to 3 percent permanent set stretch stress relief 
prior to any precipitation heat treatment. 

The soft, as quenched 0 Condition can be preserved 
by refrigeration in order to minimize springback and 
increase the general ease of forming operations. Maximum 
holding times which will preserve the formability of this 
condition are shown in Table 2 (31, Code 3206, pp. 3). 


Temperature °F 

R.T. 

32 

20 

Time 

2 hours 

2 days 

7 days + 


TABLE 2: Holding Times for 6061 in 0 Condition 

3. Specified Mechanical Properties (31, Code 3206, pp. 2) 


TABLE 3: Specified Mechanical Properties. 


Condition 

0 

T4 , T451 

T6 , T651 

Thickness - in. 

0.250 

0.250 

0.250 

Ultimate Stress in 
Tension min. ksi 

22 

30 

42 

0.2% Offset Yield 
Stress in Tension 
min - ksi 

12 

16 

35 

Ultimate Elongation 
2 inch min - % 

18 

16 - 18 

0 

H 

1 

U> 

Ultimate Shear Stress 
Typical - ksi 

12 

24 

30 
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4. Mechanical Properties ; 

In welding aluminum, it must be remembered that 
the metal in the area of the weld will be in an annealed 
condition after welding. There will be a corresponding 
loss of strength unless heat treatment is done after 
welding. Welding does not always reduce the strength 
of the heat treated alloys to that of the fully annealed 
condition because there is some air quenching of the metal 
as it cools (32, pp. 6.33). Table 4 (32, pp. 6.34) 
demonstrates this effect for the type of filler metal 
utilized in this experiment. 


Filler 

Wire 

Specified minimum 
tensile strength of 
base plate in ksi. 

Average tensile 
strength across 
weld in ksi. 

Average bend 
free elonga- 
tion % 


Not heat treated after welding 

4043 

42.0 

27.2 

16.0 

bh 

Heat treated and aged after welding 

mol 

42.0 

43.5 

11.0 


TABLE 4: Strength and Ductility of Welded Butt Joints in 

Aluminum (TIG and GMA with argon) . 


Brungraber and Nelson (33) report that 6061 alloy may 
be heated up to 550° F by butt welding without appreciably 
affecting the mechanical properties of the base plate. 

They further report that for thicknesses less than one 
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inch, the distance from the center of the weld to the 
edge of the heat affected zone (HAZ) may be predicted 
as : 

b h = 0. 6t + 0.00059 || (32) 

where E has units of volts, I amperes, V inches per 
minute, and t (thickness) inches. As an example, for 
E « 20., I ■ 240., t = .25 and V = 32.16, El/Vt - 597 
watt - min/in , and b h - .50 inch. The extent of the 
fusion zone may depend somewhat on the parameter El/Vt 
but it is governed primarily by the thickness t r.f the 
plate and the geometry of the edge preparation. 

Figures 4 through 10 provide mechanical properties 
at both room temperature and elevated temperatures. Figure 
4 provides typical tension and compression stress strain 
curves for the experimental plate. Figure 5 provides 
the average stress strain curves for varying temperatures 
from room temperature to 700°F based on a half hour 
exposure to the temperature in question. This may or may 
not be a valid approximation for the welding problem, but 
it is the only information currently available. Figure 6 
shows the effect of this temperature on the 0.2% offset 
yield stress. This may be approximated as the cubic 
polynomial given in Table 5. Figure 7 demonstrates the 
effect of the half hour soak temperature on Young's 





modulus. This curve was generated from Figure 5. Its 
cubic approximation is listed in Table 5. 

Figure 8 shows the room temperature tangent modulus 
at various stress levels. Its straight line approximation 
is 

H* x 10~ 3 (ksi) = 0.917c {ksi) + 38.064 

(33) 

The highest stress level shown on this figure is 40.7 ksi. 

4*3 

This equates to H' « 0.746 x 10 ksi and may be utilized 

effectively for constant strain hardening. Figure 9 
represents an estimate of the temperature effects on the 
tangent modulus. If one assumes that 


H' (T°F) « E (T°F) 


(34) 


then 


H' (T°F) 


H T (R.T.) x E(T°F) 
E (R.T. ) 


(35) 


Masubuchi and Andrews (11) define the constant strain 
hardening parameter f^r their one-dimensional computer 
program as 

H' " mE (36) 


From the above equations, 


m * HVE = H' (T°F)/E(T°F) = 0.0742 (37) 

This is the first such value ever determined for their 
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program. 

Figure 10 also demonstrates a computer algorithm 
for the tensile properties at various temperatures based 
on the half hour soak. To determine this approximate 
value, proceed as follows: 

a. Determine Oy(T°F) from Table 5 or Figure 6. 

b. Determine E (T°F) from Table 5 or Figure 7. 

c. Determine H' (T°F) from equation (35) and Table 5 
or Figure 9. 

d. Plot o y (T°F) at 0.2% offset. 

e. Draw line through point d above with slope 
H’ (T°F) . 

f. Connect line in e above with line through origin 
with slope E (T°F) . 

5, Physical Properties 

Figures 11 through 14 show the effect of temperature 
on the coefficient and mean coefficient of linear expansion 
a, density p, thermal conductivity k, and specific heat 
C. The coefficient of linear expansion shown in Figure 11 
is based on empirical equations for high purity aluminum 
and alloy constants provided by reference (35). The alloy 
constant for 6061 TO is 0.990, i.e., 0.990 a (pure 
aluminum) - a (alloy) . The value of this constant for 
heat treated tempers is approximately 0.015 greater. 

With these tempers, application is restricted to temper™ 
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atures which do not appreciably exceed those used in the 
final aging treatments (320-350°F). The basic alloy 
constant is limited to temperatures below 600 °F. Figure 11 
is extrapolated past 600 °F to the melting point. Figure 12 
provides an estimate of density provided by reference 
(36) extrapolated to the melting point. Figure 13 provides 
an estimate of thermal conductivity. Reference (37, 

Volume II, Part 2, pp. 769) provides thermal conductivity 
information of 6000 series aluminum alloys from 300 to 
650°K. A linear relationship was plotted parallel to 
this data passing through the one data point provided by 
reference (38). Reference (37) (Volume I, pp. 11) listed 
the specific heat- temperature relationship for pure aluminum. 
A parallel relationship (39, pp. 180) passing through the 
data point provided by reference (31, Code 3206, pp. 1) 
resulted in the estimated data shown in Figure 14. 

These physical values may be approximated by the 
cubic polynomials given in Table 5 for computer simulation. 
Appendix A lists the computer program for determining 
the coefficients given in Table 5 and for testing these 
coefficients in 10°F increments. 

Table 6 summarizes the physical and mechanical pro- 
perties of 6061 T6 in the form utilized by Masubuchi and 
Andrews (11) in their one-dimensional program. 
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0.097 
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31.70 
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5.00 
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B 
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14.26 
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16,*3 
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0.07*2 
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0,0742 

B 

0.07*2 

0.07*2 
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- 

0.07*2 

0.0742 

0.0742 


Tabla 6: Suaeeary of Mechanic* 1 and Physical Properties 

for 6051 T6 Alualmta Alloy* 
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III: PROCEDURES 


A. Scope of the Research : 

A series of six experiments measuring temperature 
and strain changes during well ;.g was performed. System 
models representing constrained bead-on-plate and butt 
joints were constructed from 6061 aluminum alloy in the 
T6 condition. Welding procedures were utilized which 
offered full penetration with minimum heat input. This 
allowed strain gage locations as close to the weld 
line as possible without exceeding the 400°F maximum 
allowable gage temperature limitation. 

B. Selection o f Parameters : 

The aluminum selected for this investigation was 
picked because of its easy availability and wide use in 
the marine and aeronautical industries. It develops its 
strength from its heat treatment. The 0.250 inch plate 
thickness used in this study was determined by availability, 
past experimental work, ease of handling, and an interest 
in approaching a two-dimensional problem, yet minimizing 
the extent of bending and buckling. The joint design 
for the butt welds used in this experiment (straight) 
was selected for simplicity and ease of preparation. 

Steel backing plates were utilized below the weld and at 



the clamped ends of the plates to provide support and 
distributed clamping and yet minimize heat conduction 
effects. Each specimen was clamped at its edges to a 
0.500 inch mild steel bed plate which resisted deformation. 

The size of the test plates was set at a nominal 18 
by 30 inches to provide essentially steady state conditions 
at strain gage locations, and to approach the temperature 
distribution of an infinite plate. 

The weld process used was semi-automatic Gas Metal-Arc 
(GMA or MIG). It allows control of weld variables, reduces 
operator error, and fosters repeatability. Ho preheating 
was necessary. 

C. Strain Measurement by Electric Resistance Strain Gages: 

The fundamental concept of strain gage operation 
is that certain conductors exhibit a change in electrical 
resistance with a change in strain. Gages designed 
according to this principle are attached to test materials 
whose strains are then monitored by measuring resistance 
variations across the gage. In the case of welding 
thermal strains, the observed resistance change, is 

made up of: 
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where 


AR l l£ e> 


A W 


AR 3 (aT) 


ar 4 (T) 


the resistance change corresponding to 

elastic mechanical strain, e , from which 

e 

stresses can be computed* 

the resistance change corresponding to 

plastic mechanical strain, e , if it exists. 

hr 

the resistance change corresponding to 
temperature induced thermal strain, aT. 
the resistance change caused by thermo- 
electric effects in the gage itself. 


While it is not presently possible to discriminate 
between the two mechanical strains, z^ and Ep, AR^faT) 
and AR 4 (T) can be separated out by physical temperature 
compensation of the strain gage material coupled with 
empirical calculation. For this investigation, a T for the 
gage was approximately 13.0 microstrain/°F. The AR 4 (T) 
and the difference between gage a T and base plate 
will then generate an error. To compute this error, 
a test gage of the type and lot used in the experiment 
is mounted on a small sample of 2024 aluminum alloy by 
the gage manufacturer. The sample is then heated at 
equilibrium until a curve of "apparent strain'* vs. temper- 
ature is obtained to the operating temperature range. 

The gage readings recorded in the weld experiments can 
then be corrected by subtracting out the apparent strain 


... 1 v. rxxxasSLs. i r *Z ' 


value corresponding to temperatures observed at the gage 
location. Figures 16 and 17 show this factory provided 
information for the individual gages and rosettes 
respectively, in each case, the curves have been 
reduced by the factory to fourth order polynomials yield- 
ing microstrain for ease in computer calculation. 

Gage: 

EAP = -*69.36 + 2.07T - 2.02 x 10~ 2 T 2 + 5.39 x 10" 5 T 3 

- 3.27 x 10“ 8 T 4 (38) 


Rosettes: 

EAP - -60.68 + 2.22T - 2.40 x 10“ 2 T 2 + 7.05 x 10~ 5 T 3 


- 5.59 x 10“ 8 T 4 


(39) 

Since 


6061 T6 has a greater coefficient of thermal expansion 
than the 2024 T4 test plate, a correction to these poly- 
nomials has been made. This correction was calculated 
as follows: 


In general, 

Thermal Strain - 

From reference (35) , 


a (T) dT 


(40) 


a 6061 ~ C 6061 a pUre A1 * = ,9S a pUre Al * * 41 * 

= .97 a pure Al. (42) 

(43) 


a 2024 C 2024 a pure A1 


a pure Al. = 
T* 


6061 

.99 


■ !. 


a 


pure Al . 


AEAP 


6061 ” c 2024> 


(T) dT 


(44) 
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*2 

AEAP “ | <f§£> ta 6061 (T)) dT 

T 1 

Integrating the polynomial given in Table 5 for 
0 £ T £ 400 Q F, between 77 °F and T°F yields the following 
apparent strain correction for both gages and rosettes; 

AEAP = -15.4979 ■+ .2463T + .9561 x 10 _4 T 2 - .7879 x 10~ 7 T 3 
+ .5303 x lO -10 ^ (46) 

Equation (46) must be added to equations (38) and (39) 
to obtain true apparent strain corrections. The true 
apparent strain is then subtracted from the actual strain 
to provide true mechanical strains. 

Strain gages manufactured in the form of rosettes 
provide as many as three independent readings in three 
directions at a single location. This is enough to 
completely describe the two-dimensional strain state 
at that location (41, 42, 43) and is discussed in the 
next chapter. 

D. Description of Apparatus; 

1. Test Plates: A description of the test 

specimens is summarized in Figures 18 through 22 which 
provide dimensions and gage locations. Devices^) through 
(J^jwero dynamically measured during welding. Gages 
(5l) through were statically measured before and after 




Figure 19: Test No* 3 (Bead -on- Plate Weld) Measuring Temperature 

Distribution and Strains. 
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Figure 21: Teat No* 5 (Butt Weld) Measuring Temperature 

Distribution, Strains, and Extension of Root Gap* 
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Figure 22: Test No. 6 (Butt Weld) Measuring Teaporcture 

Distribution, Strains, and Extension of Root Gap. 
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welding. The arrangement of constraining clamps is 
shown in Figure 23. 

2. Sensors and Instrumentation ; 
a. Strain Gages. Two types of gages were used 
in this investigation: SR-4 foil 45° rosettes and 

SR-4 foil single element gages. Gage properties 
were as follows: 


Gage 

Designation 
Manufacturer 
Grid Length (in.) 
Grid Width (in.) 
Overall Length (in.) 
Overall Width (in.) 
Temperature Range 
Resistance 
Gage Factor 
Cement 

Protective Covering 


SR-4 Rosette 
FAER-18PB-12S13-ET 
BLH Electronics 
3/16 
0.90 
0.280 
0.540 

-50 to +400°F 

120 . 

2.03 
EPY-600 
BLH Barrier-C 


SR-4 Gage 
FAE-25-12S13L 
BLH Electronics 
1/4 
0.13 
0.35 
0.13 

-50 to +400 °F 
120 . 

2.07 
EPY-600 
BLH Barrier-C 


b. Temperature Sensors: All temperature sensors 

used were Chrome 1/Alumel thermocouples made from 
Leads and Northrup No. 28 wires. Each thermocouple 
was spot-welded onto the test specimen and protected 
by No. 33 Sauereisen Sealing Cement.' Each was 










located 0.536 inch (1 second of arc travel) ahead 
of the strain gage mid-axis. 

c. Instrumentation. Strain gageii were connected 
into a Potentiometric Circuit (Half -Wheatstone 
Bridge) , balanced and calibrated as indicated 
schematically in Figure 24. Thermocouples were 
referenced to a 32°F ice-bath and calibrated as 
indicated in Figure 25. Both circuits were fed into 
a Honeywell continuous-recording, 12-channel 
Visicorder. When the raw data was actually read 
off the recorder tape, some traces were delayed or 
advanced with respect to others. This was done 

to correct for the finite difference in position 
along the weld line of the thermocouples (1 second) 
and strain rosette elements (+ 1/3 second) . 

ThiG produces the effect of a simultaneous 
strain and temperature reading at the strain gage 
location. Timing was accomplished by means of an 
electric stop watch as well as a timer integral 
to the Visicorder. Static strain measurements 
were taken with a Budd Instrument Division Model 
P-350 strain indicator. 

d. Extensiometer. Figure 26 shows a sketch of 

the extensiometer utilized in Tests 5 and 6. Strain 
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Battery 


Figure 24s Strain Gage Inatrunentatlon Circuit (7.) 
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Flgura 23: Tharaocdupia Instnmantatlon Circuit i 7). 
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gages were located on both sides of a spring steel 
ribbon and connected into adjacent legs of the 
resistance network circuit described in paragraph 
c. above. The net effect of such a set-up is to 
double the bending sensitivity and null the stretching 
sensitivity. Figure 27 provides the calibration 
curve for the Visacorder referenced to actual 
dimensional changes in inches. The extensiometer and 
its mounting tabs may be seen in the photographs 
in Figures 28 through 31, It was pulled back away 
from the torch to prevent damage when the arc past 
the mounting tabs. Once the torch was safely clear, 
it was set back in place to observe weld shrinkage. 

It is noted that the mounting tabs were cemented in 
place with EPY 600 cement vice bolting or welding 
to minimize local straining of the base plate. 

3. Welding Equipment and Conditions ; Welding conditions 
are summarized in Tables 7 and 8. The welding machine 
utilized was manufactured by the Linde Division of Union 
Carbide Corporation and consisted of a type HW-16 GMA 
torch, a SVI-300 power supply, and associated governor, 
carriage and wire feed mechanisms. Travel speed, arc 
voltage, amperage and arc length were present to the same 
value for all experiments. Some fluctuations of amperage 
did exist when the geometry of the weld changed .from 




f | RU re 2«: Overview of Exper 1 rrent* 1 Equipment Showing Tnetruirentet « on end Recorder. 






Figure 29: 


Plate Instrumented and Clamped Under Welding Torch. 
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H ?ure 30: Piste Instrumented .ind ripped LY.der Welding Torch 





bead-on-plate to butt. Wire feed is a functin of the 
other variables and was maintained by the machine auto- 
matically* Figures 29 through 31 are photographs of 
the actual equipment set-up. 

E. Experimental Procedures: 

The experimental operation is shown schematically in 
Figure 32. The test plate was instrumented and clamped 
into place. The welding machine was lined up with the 
joint and positioned over the run-off tab at tho left 
end of the backing plate. Welding speed, arc voltage and 
amperage were pre-set. The visicorder was actuated and 
an arc was struck on the backing plate. As the welding 
torch crossed from starting tab to plate, the timer was 
started. The Visicorder output was marked when the arc 
passed the center strain gage location. When the welding 
bead reached the run-off tab at the right end of the 
plate, the arc was extinguished and the plate allowed to 
cool. Thr recorder continued to monitor the gages until 
the plate cooled to ambient temperature. After cooling, 
the plate was released from its clamp, and the new strains 
recorded either with the recorder or the Budd instrument ^ 




•Figure 32: Schematic of Apparatus and Procoduroc (7) 


IV: DATA DEDUCTION 


A. G e neral: 

At any point on a free (unloaded) surface of a 
solid, it is necessary to know three independent quantities 
in order to specify the state of stress completely. These 
quantities are the magnitudes of two principal stresses, 
a^, and o 2 , and their directions, <J> or (<f> + 90°), with 
respect to some reference. For isotropic elastic materials, 
these values can be calculated from strains measured on 
the surface at the point in question, and since three 
independent quantities are to be determined, it will be 
necessary to make three independent measurements of strain 
(42, pp. G2-03) . 

B. Strain Calculations; 

Equation 47 describes the relationship for general 
strains for a three- element rosette shown in Figure 33. 



Figure 33 
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In the special case of a three-element 45° rectangular 

rosette shown in Figure 34, 
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Equation (47) reduces to 
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Solving Equation (48) for e . e , and y , 

x Y xy 



(48) 


(49) 


The Mohr's Circle for this special case is shown in Figure 35. 
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Figure 35: Mohr*o Circle for the Re<tangular Rocette With 

Three Observations of Strain <42, pp. G2-12). 


Principal strains arc 
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e ) 2 + 

y 



(50) 


E 


2 





V 2 + 



(51) 


where ( f ) is defined as the angle of one of the principal 
axes with respect to the axis of reference. Mathematically, 


Tan 2cj> 



(52) 


In the literature, there is considerable confusion 
concerning t J> and which principal axis it defines. 

Murray (42,. pp. G2-11, 15) offers a consistent method 
for dealing with this problem. If 4»^is defined as the 
angle measured (positive in the anti-clockwise direction) 
from the positive OA axis of the strain rosette to the 
positive 01 axis which corresponds to the direction of 
e^, the algebraically largest principal strain, then: 

e + e 

Rule 1. When e> — , then 0° < 4> < 90° (53) 

b 2 i 

e + e 

Rule 2. When e < — - , then -90“ < <f> < 0° (54) 

b 2 

e + e 

Rule 3* When e. = — - ( 55 ) 

b 2 

and (a) £ > c then e = e and = 0° (55a) 

a c a x J- 

or (b) e < e then e ^ e and 4> . = + 90° (55b) 

Q. C Q u J- 
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Rules 1 through 3 are defined with axis OA as the 
reference. To convert to the situation consistent with 
this investigation, axis OA is t 45° (anti-clockwise) from 
the x-axis. 


B 



* x " ^1 + 45 ° (56) 

The geometry utilized in this investigation reduces 
to that of the mathematical model if the axes are kept 
consistent with Figure 36. 

C. Stress Calculations : 

This investigation yields mechanical strains. The 
automatic temperature compensating feature of the strain 
gage coupled with the apparent strain correction removes 
the thermal strain effect. Stresses must be determined 
by total strain. 

E mech + ^thermal E total (57) 
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where T 


G the 

rmal j T 

a (T) dT 

1 


(40) 

For 

general two 

-dimensional 

stress calculations 

9 

a 

X 

E 

1 - V 2 

( e + 

x total 

VG ^ total’ . 

(5B) 

a 

y 

E 

1 - v 2 

“W + 

ve x } 

-total 

(59) 

c i 

E 

1 - V 2 

(e. + 

total 

) 

total 

(60) 

°2 

E 

2 

1 - ^ 

(e 9 + 

total 

ve, ) 

’’‘total 

(61) 

For this 

particular 

configuration, 


v « e - 

x ^mech a 

• e = y 

c ^total 

(62) 

T 

xy 

= Gy - 

’xy 

E 

2(1 + v')" 

Y X y 

(63) 


v is Poisson's ratio. Nominally, ^ K .333 for aluminum. 

^ will probably change with temperature, but should not 
exceed a maximum value of .500 in the weld puddle because 
of volume consideration when melted (at .500, there will 
be no volume change) . 

Stress calculations are dependent upon E(t) and a(T) 
assuming v is relatively independent of temperature. These 


T 


T 


1 
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values may be determined from Table 5 in Chapter II. 

The values of E(t) given assumes that the value from a 
half-hour soak approximates that from welding. 


D. Determination of Plastic Conditions: 

Masubuchi (43) theorizes that the region in the 
vicinity of the weld has undergone plastic deformation 
when the invarient I is larger than the value of yield 
stress/ where, 

I = {a 2 - 'j a + • a 2 + 3t 2 ) X ^ 2 ( 64 ) 

x x y y xy ' ' ' 

This value was derived from the second stress deviator 
tensor invariant (44, pp. 41) which for a two-dimensional 
situation is reduced to 



a a + cr 2 + 

x y y 



§ i 2 (65) 


The Distortion Energy Theory (Von Mises r Yield Criterion) 
at the yield point in simple tension (44, pp. 75) predicts 


that 



( 66 ) 


Equating Equations (65) and (66) 

12 1—2 

J 2 = y 7 = fo < 67 > 

or 

1 = 0= /iJj (68) 
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Therefore, a check to see if plastic conditions exist 
would occur is that 

l-Q) > 1.0 (69) 

a *T) 

It is important to note that any single component of the 
stress field in Equation (65) may be greater than the 
yield stress in simple tension without having yield occur. 

Again, the yield stress is determined from Table 5 
in Chapter II. 


Data Reduction Program: 


Appendix B lists the basic program utilised for 
rosettes in this investigation. It is modified slightly 
to suit other geometries (rosettes in different locations) 
or slightly different welding conditions. Input cards 
are listed as comment cards in the main program and 
samples are shown in the appendix. In the case where all 
strains are zero or equal, the program defaults to $ - 90°. 
For example, if at time =0.0 the strains are initialized 
to zero, the resulting 4> at time = 0.0 is meaningless. In 
the case where initial strains at time = 0.0 are not zero, 
the program will reference all succeeding strains to 
their initial v. oie. Single element strain gage data 
may be processed by this same program to provide temperature 
compensation. 



V RESULTS : 


A. General Trends: 

Figures 37 through 59 graphically present part of 
the experimental data tabulated in Appendix C. Table 7 
discusses the particular problems and parameters associated 
with each test. Several effects are discussed below: 

1. Arc Efficiency . Figures 37 through 40 show the 
temperature distributions for Test No. 1 and 2, and also 
the one dimensional computer analysis results for 0 = 0.70 
and 0.75. The n = 0.75 curve seems to fit the data best. 
This efficiency factor includes all heat loss effects 
associated with the heat generation and cooling. The 

tail of the experimental temperature distribution is higher 
than the computer analysis because , geometrically, the 
experimental plates are not infinite and cool slower than 
the computer model. 

2. Placement of Strain Gages . Figure 41 shows the 
transverse temperature distribution based on Figures 37 
through 40. In general, the peak temperature predicted 

at transverse locations are supported by the testing shown 
in Figures 42 through 44. Figure 41 provides an effective 
technique for positioning strain gages to prevent damage 
or exceeding their temperature compensation limits. 


Test 

Number 


I 


V 


Amperes 


Volts 


V 

Inch /Sec. 


Heat Input 
Joules/Inch 


Comment 


1 . 


230 


20 


0.536 


8582 


Butt weld to measure temperature distribution only. 
Immediately opposite thermocouple line arc burned 
a deep crater in the weld. Temperature distribution 
at 0.675 inch reflects problem. During last quarter 
of run, arc cut more plate than it welded. More tacks 
are required, i.e. , shift from 2 to 4. 


2 . 


3. 


4. 


220 

240 


240 


5. 


6 . 


250 


250 


20 

19 


19 


19 


19 


0.536 

0.536 


8209 

8507 


0.536 


8507 


0,536 


0.536 


8862 


8862 


Bead-on-plate to measure temperature distribution only 

Bead-on-plate to measure temperature distribution 
and strains. Sensitivity of Visicorder set low. 

Bead-on- plate to measure temperature distribution 
and strains. Strain readings in the vicinity of 
the arc on all rosettes appeared erratic. Strains 
were much higher than expected, 

Butt weld to measure temperature distribution and 
strains. 

Butt weld to measure temperature distribution and 
strains. Gages 4 and 5 failed after arc passed 
by in vicinity of highest strains at 30+ seconds. 

At 42 seconds, Visicorder paper ran out and all 
remaining data was lost. 


Table 7: Summary of Welding Conditions for Each Test 


C-i 

.ts. 



























Test Number 

I: Amperes 

V: Volts 

v: in. /min 

t: inches 

El/vt: Watt-min/in^ 

b^inches 

1 . 

230 

20 

32.16 

.250 

572.14 

0.488 

2. 

220 

20 

32.16 

.250 

547.26 

0.473 

3 * 

240 

19 

32.16 

.250 

567.16 

0.485 

4. 

240 

19 

32.16 

_ 

.250 

567.16 

0.485 

5. 

250 

19 

32,16 

.250 

590.80 

0.499 

6. 

250 

19 

32.16 

. 250 

590.80 

0.499 


TABLE 8: Summary of Extent of HAZ (b^) 
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3. Butt vs . Bead-on-plate Characteristics . Figures 41 
through 44 provide temperature distribution data for the 
tests with strain gages installed and Figures 45 through 
51 provide data on transverse and longitudinal strain dis- 
tributions resulting from these tests. Ideally, a butt 
weld becomes more and more like a bead-on-plate weld for 
multipass situations. One would hope for a two dimensional 
model which cculd cover both cases. Unfortunately, 

Figures 43 and 44 show a different inicial temperature 
distribution for bead-on-plate and butt welds. This implies 
that for the plate thickness utilized in this study, the 
two types of weld will behave somewhat differently in 
the way of thermally-induced strains, moments, and distor- 
tions. The bead-on-plate weld is not a two dimensional 
situation, but the butt weld clearly is. Figures 45, 46, 

48, and 49 seem to agree with this observation for strains, 
since bottom bead-on-plate values do not match the top ones. 
The higher transverse strains for Test No. 4 are noted, 
but not fully explained. The strains were eratic at their 
peak values as if the plate was slipping in its constraint. 

4 . Correlation Between Experiments . In general , 
there is not much correlation between residual strains from 
experiment to experiment. This problem has been noted at 
MIT with previous experiments and is again noted here in 
the various strain distributions. An effort to find a 


means of correlation was made by looking at the angle 
between the largest principal strain and the x-axis (4> ) . 

A 

This data is shown in Figures 52 through 55. There seems 
to be little correlation in the residual strain region 
other than that the major strains are either close to 
the x or the y axis. Along the edge of the panel, however, 
the test data consistently show a <j> which implies loading 
by shear. 

Figures 56 through 58 provide a much more successful 
means of correlation. The ratio of the invariant described 
in Equation (64) to the temperature dependent uniaxial 
yield stress from Table 5 is plotted with respect to time. 
For residual conditions, there appears to be very good 
correlation between the "state of stress^ composed of 
all stress components from test to test. In addition, when 
the value of this ratio exceeds 1.0, a yielding situation 
exists, and plastic deformation occurs. 

The static residual strain measurements made during 
this study offered such a wide and inconsistent variety of 
values that they were not included. One possible reason 
for this is that the position of the strain gages one inch 
away from the centerline of the plate may have been near 
the region in Figure l - 4, where the stress distribution went 
through a zero point with a large slope . Hence, a minor 
variation in stress could imply large positive or negative 
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variations in strain. 

5. Extensiometcr . The extensiometer results are 
plotted in Figure 59. Initial measurements vary from 
experiment to experiment because of the difficulty of 
placing the holding tabs in the exact same position. 

However, it is interesting to note that for the two 
situations presented here, once the weld was filled in 
the vicinity of the extensiometer, the extensions became 
a measurement of transverse shrinkage and were similar in 
both tests. Test No. 1 utilized two tack welds at the 
ends of the plate, and the arc became unstable as a result 
of metal movement in the center. All subsequent tests 
utilized four equally spaced tacks, so the plates were 
essentially rigid except in the immediate region of the 
arc. 

B. Accuracy of the Experimental Model and Instrumentation : 

1. The Physical Model . As mentioned previously, 
the plate geometry, level of constraint, and welding pro- 
cedures are considered to be representative of ship 
fabrication processes. However, the means of constraining 
the panels is unsatisfactory because of evidence of slipping. 
Considerable difficulty was experienced in cementing the 
strain gages to the model. The best cement still remains 


EPY-6QQ because it offers a lower curing temperature. 
However, pressure must be applied to the gage while the 
cement is curing, making this an extremely critical stage 
in model production. 

2. Data Reduction Calculations . The temperature 
compensation has been previously discussed. However, in 
the strain distributions presented, there is a strange 
fluctuation in the results in the vicinity of the arc 
immediately before the thermal strains take off. It is 
of the opinion of this author that this is in part, if 
not all, due to temperature compensation techniques and 
represents only an "apparent strain". In any case, it is 
a minor effect and relatively unimportant. 

3. Instrumentation * The Visacorder output was 
scaled initially in Test No. 4 to 1200 microstrain/inch 
and for later tests at 500 to 600 microstrain/inch. The 
reading accuracy is well within the 5% instrument accuracy 
estimated by technicians aiding in this study. 

4. Welding Machine . The welding machine utilized 
in this study leaves much to be desired. The speed drive 
appears not to give a uniform arc velocity although the 
time to span the full plate length seems to be consistent. 
The voltmeter and ammeter on the power supply have never 
been calibrated and a clamp-on ammeter and separate 



voltmeter had to be utilized as a result. The meters on 
the actual welding machine controls are either inoperative 
or way out of calibration. There is no screen over the 
tip of the torch which prevents back splatter from fouling 
the automatic filler metal feed system. Luckily, the wire 
jams which occurred as a result were on the starting tabs 
and the expensive test plates were not ruined. 
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VI RECOMMENDATIONS : 


1. The small number of tests performed in this study 
should be expanded into more tests, ultimately with plates 
of varying thicknesses. Data reduction from the Visa- 
corder strip charts is difficult and time-consuming. If 

a large number of test are run, it becomes almost imperative 
that some form of direct digital output be utilized to 
provide punched tapes or cards compatible with MIT's IBM 
computer system. Multipass tests likewise require this 
capability. Voiding is a statistical process, so the 
larger the basis of data the better the results will be. 
Obviously, there are serious questions of financial and 
manpower resources available, so the trade-off of more 
experiments vs. better data must be seriously considered. 

2. More bead-on-plate tests which record top and bottom 
strain and temperature distributions are needed to verify 
a two dimensional approximation. 

3. The welding machine utilized by tV • study should be 
carefully inspected and the discrepancies noted in Chapter 
V remedied, or another machine must be found. 

4. A "V" or "J" weld preparation should be considered 
for butt welding thicker plates. 
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5. Smaller strain gages should be utilized to allow 
easier application and closer positioning to the weld. 

This sacrifices the larger area for averaging measurements, 
however. Multiple-element gages should be of the stacked 
variety to eliminate time corrections for individual 
elements, and make data reduction simpler. Gages such 

as BLH Electronics' FABR-12-12SX13ET (stacked, three-element 
45°, 0.12 inch gage length), FABX-12-12S13ET (stacked, two- 
element, 90°, 0.12 inch gage length), or FAER-06RB-S13ET 
(three-element, 45°, .0625 inch gage length) or equivalent 
Micro-Measurements, Inc. gages should be utilized. 

6. If time permits, an independent check of apparent 
strain corrections is recommended to verify factory data. 
This may be accomplished by cementing a gage to a sample 
test plate and heating it in an oven. 

7. High temperature strain gages should again be con- 
sidered to see if they are any more reliable th«n reported 
previously by Klein (7) . 

8. The test plate mounting assembly should be modified 
to allow bolting vice clamping of e plate boundaries. 

This may be done simply by drilling and tapping the bed 
plate currently in use, and will eliminate any slipping 
from taking place. 
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9. Further use of the oxtensiometer is not recommended. 

It costs one channel on the Visacorder and causes the 
loss of a rosette (12 channels are available: 3 rosettes 

and 3 thermocouples, or 2 rosettes, 2 single-element gages, 
the extensiometer and 3 thermocouples) . This author feels 
that there is more to be learned at present from the "state 
of stress" than from shrinkage measurements. However, if 
continued use of the ex tensiometer is desired, two more 
gages should be installed to double its sensitivity. 

10. Continued basic research is needed to provide better 
physical and mechanical properties for a larger variety of 
structural materials. 
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C MICRCSTRAIN. 

RFAL m 
IN TFCFR GAGE 

DIMENSION T(99> ,S2ER0(99, 3), $TRN(99*3J ,SC99, 3) 

DTMFNSICN EXI991 ,EYf99) ,GXYI99 ) ,El (99) ,F2( 99 ) 

DIMENSION PHI X (°9 ) 

HI MENS I CN 6(4) ,E(4)?Yt 4>,ETH0(?)* SIGXt 3 » , S T GY{ 3) ,T AUXY( 21 • X4 ( 3 ) 
DTMFNSICN Y1I3) ,X*3f 3) 

READ (5, SO) GAGE, VOLTS, AMPS 
C C FOR w AT ( 15, 7*6.9) 

JsGArr 

WRIT Ft 6, 100) 

t FORVATCl^X, ‘WELDING STRAIN ROSETTE CAT A: EXPERIMENTAL RFfUlTS*/ 


T5X, ************ *********************************** **** ************ 
2*'///) 

WRITF(6«l0l ) 

101 FORMAT! 18X, * TEST NUMBFP 4: BE AD- CN- PL ATE*/) 

WR I T E (6 * 110 1 VOL T $ , A KPS 

110 FORMATC 18X» ’ALUMINUM ALLOY 6061 • IN T65 1 CONDIT ION* / IflX *• THICKNESS 
1* 0.250 INCHFS ' /1 8X * •GMA WFLD PROCESS* /18X.» ARC VOLTAGE - »,F6.0/ 

7 1 8 X. * AMPERES = * «F6.0/IBX*$PEEQ CE TRAVEL * 32.16 1NCHF S/^INLTE* ) 

W R I T E 16*111) 

111 F0RMATU8X. *4043 FTLLEP METAL: 0.C625 INCH D! AM{=T6 R* / 1 8 > , 

?■ UNITS OF STRAIN - K1CRCSTRAIN*// ) 

WRITE (6*11 2 ) 

II? E0RMAT<S3X, *CAUTI0N*,//l8Xt »TEKP = 0.00 COMPUTATIONS REPRESENT*/ 
21.8 X ♦ 1 BLANK SPACFS ON INPUT CARCS AND ARE’/tSX, 

T’MEANINGIESS. IGNORE ALL SUCH COMPUTATIONS. *///! 

WRITE! 6f!20) 

120 FORMAT! ?8X,» TRANSVERSE DISTANCE FRC* WELD IN INCHES*// 

210X » *TI ME , *18X« * i. 00* * 12X,* t.OC* *12X* # 1.G0*/ 

3RX. * I SECONDS M .11X, • (CENTERLINE I * ,4X, * ( EGTTOM C/L * * t 
46X.MEDGEI *////» 

RFACC5.SUA 
REACTS* 91 ) E 
READ t 5* 91 ) Y 
91 FORM AT ( 4FI 5.71 
READ (5* 52 )NU 
9? FORM AT(F10.3| 

READ IN ROSETTE DATA: TIME* TEMP, STRAINS A,p AND C. 

I R e A n ( 5 * 1 4 0 , END* 999 I T ! ME . C T ( K > * K * I . J 1 * U S ( K , u * L *1 * 3 1 * K= 1 * J > 

14C EQRKAT(F3.0*12F6.2) 

IFITIME.GT.O. IGF TO 3 

on f il*!* j 

ETHfni )*Ail1*T(IU*A(2|*T( IIi**?/?.4Ar?)*T( II) **V 3.* A! 4 ) *TI t T ) 
l **4/4. 

00 7 KK=l,3 

7 S7FQCU T.*K)*$| TI*KK| 

6 CONTINUE 
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3 00 8 JJ*1.J 

C TEMPERATURE COMPENSATION FOR FAER— ISR6— 12S11— FT LOT NO* 8263 

C GAGFS ON 20?4-T4 ALUMINUM TEST FL ATP* 

ASTRa-60* 6R+2*22*T f JJ )— 2*40E~2*T ( JJ I **2+7* 05E—5*TI JJ 1**3 
1-5.5<5F-F*T<JJ)**4 

ASTRl^l9*4979+.2463*TCJJl*.9S6LF-4*Tf JJ|**?-.7979E-7*TC 
U.5303E-iO*TUJ >**4 
ASTRsASTR+ASTRl 
TFCTIMF.EO.O.t GO TO U 
00 <3 KK*t * 3 

9 STRNl JJ*KK| * SC JJ*KK)-rZERO(JJ*KK)-ASTR 

GO TC 13 

1 1 00 12 KK*l « 3 

12 StRNt JJ*KK| -SC JJtKK l-ASTR 

11 CONTINUE 

C STRAINS ARE NOW TEMPFRATURF COMPENSATED * 

CALL STRAIN! STRNt JJ * IT ,STRN(JJ.2l ,STRN| JJ,3> *FX( JJI ,FY f JJ) *GXY IJ%! ) 
2* FI CJJl «E2 ( JJ ) *PHIXCJJ1 ) 

CALL STPESSCeXC JJWEYI JJ)*GXY( JJUETHO C JJ> , A, E.Y,TUJI *NU 
ISIGXCJJ I.SIGYC J.IKTAUXYCJJ J.X4CJJ) *Yll JJI.X5C JJH 
3 CONTINUE 

WRlTFC6*l60ITtMF»lTCK)tK*l *J) * ( ( SC K*L) t.K*l *J|*i =1*31*1 <STRN!K*U*K 
2*l#J)*L=i*3)*(EXCKI *K=l*41 *(FYCKJ*K=l*J 1* CGXYCK) ,K=t*J»* CFICK),K*l 
3*J) * (E2CK) *K = l* Jl* CPHIXCK) *K=l*J) 

160 FQRMAT16X* F8* 2/5X* * TEMPERATURE (DEG*FAHR J * * 3{ F10* 2* 6X1/ 

?5X*‘ STRAIN A. f»EAS»‘*5X*3|F!0. 2*6X1/ 

3SX* ‘STRAIN 8. ( MEAS 1 * * 5X*3 (F10.2,6X) / 

4SX* * STRAIN C. CMEAS»*t5X,3{F10.2,6X)/ 

5SX. ‘STRAIN A. tMECHl‘*5X*3fP10.2*6X»/ 

6*5X, ‘STRAIN P. (MFGMJ* *SX*3tFlO#2*6X) / 

75X • • STR AIN C* I NFCH 1 * * 5X* 3 IFlO .2*6X1 / 

3SX* ‘STRAIN X* f MECHI* *5X*3(Fi0.2.6X)/ 

9SX* * STRAIN Y. ! MECH M *SX*7( FlO. 2 * 6 X ) / 

15X * ‘GAMMA XY. ( MpCH » * t5X* 1 ( c 10 . 2. 6X) / 

1SX* ‘CRIN.STP.1. (" 1 • *SX*3{F10.?*6X) / 

35X* ««RT\.STR.2. C« 1 * * 5. X * 3' f F 10 . 2 * £ x ) / 
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c\r» n o n n n.n n o o o 


45 X * * PHI X (HEr.PPESI •*5X*3IF10,?.6X)) 

WRlTEff* 161 ) (STrXtKj*K-|* J|, (SIGY(K),K*1, J), (TAUXYCKJ ,K*1* J>* 
i!Yl(KUK=l, J),tX4(K JtK=l,J) ,CX5(K) ,K=1»J) 

If 1 FOP«AT(5X** SIGMA XX { KS II • * ?X , 3 ( F l C.2* 6X ) / 

15X* * SIGMA YY IK SI J * ,7X *3 f F10.2 *6X ) / 

?5X* * T AU XY (KSn* »7X,3CF10. 2*6X1/ 

35X,*YIFtn STRESS (K St > • ,3X * 3 f F 10 • 2 *6X I / 

45X, *1 IKS 1 f * *14X*3( E10.4*6X>/ 

55X* • I/YIFLO STRFSS • *7X , 3 ( E 10.4* 6X I /// ) 

GO TC l 

9<39 CALL FX IT 
ENO .* 

5U8R0UT INF STR ATN ( E A * EB , EC. * FX* FY * GXY * F I * F2 *PM I X ) 

THIS SUBROUTINE CALCULATES STRAINS ALONG THE X* Y AND PRINCIPAL 
AXES: THE SHEAR ANGLE GAMMA XY; AND PHIX* THE ANGLE BETWEEN THE 
X AXIS AND THE AXIS OF THE LARGFST PRINCIPAL STRAIN* 

CALCULATIONS AR F BASEC CN THE RECTANGULAR ROSETTE WITH THREE 
OBSERVATIONS OF STRAIN. 

PH!1 IS DEFINED AS THE ANGLE MEASURED POSITIVE. IN THE ANTI- 
CLOCKWISE DIRECTION FROM THE POSITIVE CA AXIS OF THE STRAIN 
ROSETTE TO THE POSITIVE 01 AXIS WHICH CORRESPONDS TO THE DIRECTION 
OF THE LARGEST PRINCIPAL STRAIN. 

COMPUTATIONS* OTHER THAN PHI, UTILIZE AXIS OX AS THE REFERENCE 
AXIS, 

THE ANGLE BETWEEN THE X AXTS ANC THE OA AXIS = 45 DEGREES. 

ex=ea-eb>ec 

FY-EP 
GXY^EA-EC 
A=EX4EY 
B« FX-EY 

C-SCPT(R**?+GXY**2) 

A+C)' 

E «-C» 

Bt« c A-EC 
B2= EA+EC 
B3»2‘**FR-B7 


i 
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IF LSI. SC. 0.) GO TO fc . 

' PHI l*A8S( .5*ATANIB3/Bl II 
GO TP 7 

6 PHli *3.1416/4. 

7 : CONTINUE * 

D*,5*tEA+EC* 

IF I FP.GT.O JGO TO 1 * 

Iff FE.LT.D |G0 TC 2 
IF IE A ,GT. EC I GO TO 3 
IFIEA.LT.EC1GO TO 4 
l PHI2*+PHII 

GO TO 5 

7 PHI2*-PHII 

GO TC 5 

3 PH I 2*0.0 

GO TC 5 — 

4 PHI2*+3.1416/?. x 

5 PHI 3*PH I?*| 80./ 3. 14 16 
PHIX=PH 13+45.0 

Q 45 CEGRFE CORRECTION IS MADE TO BRING THF REFFRENCE AXIS 
C FROM PA TO X* 

RETURN 
END \ . 

* SUBROUTINE STRESS ( FX, EY,GXY,FTHO, A ,F, Y,T ,NU ,S!GX» SIOY,TAUXY. 

’ mm . ra 1X4,VL*X51 

I:?!. real nu 

DIMENSION A (4) « Ef 4 ) , Y( 4 ) 

[••: C AT THIS POINT, THERMAL STRAINS ARF CALCULATED ASSUMING ISOTROPIC 

C MAT FPIAL. 

FTRT*Af l ) *T+AI2 J*T *♦?/?•+ A | 3 J *T**3/3. + A 14 )*T**4/4. 

ftf*etht-etho 

r AT THIS POINT, TOTAL STRAINS ARE CALCULATED, 

rxT*EX+FTH 
EYT*EY*ETH 

Fl-Ffl V+FI 2 >"*T + E|3 |*T++?+F| 4l*T*+3 
r Ft yniNG+S MODULUS at TEMPER A TijRE T, 



9ZT 
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Intentionally Left Blank 



1 




Xt=1.0F-3*El/U.~NU**2> 

X?-1.0E-3*El/(2*n .*NU) » 

SIGX-X1 *t EXT-t-NU A EYT I 
$IGY*Xt *1 EYT*NU*€XTi 
TAUXY=X?*GXY 

C UNITS OF STRTSS ARE NOW IN KSI. YOUNG* S MnOUlUS WAS READ IN 
C AS PS! * F— 6 AND STRAINS WERE IN TFR«S OF MICRO STRAIN* XI AND 

C X2 HAVE BEEN CORRECTED BY E-3 TO CONVERT PST TO KSI. 

X3 = <Sfr.X**2-SIGX*SIGY+S!GY**?+3.+TAUXYI 

X3 l* APS (X3 I 

XA*S£RTtX3ll 

c K4 is HASUPUrHI'S STRESS INVARIANT IN LNTTS OF KSI. 

Yi^yaiTYU i*t+vm*T**2*Yt a>*t**i 
t Yl IS YIFLO STRESS AT TEMPERATURE T IN UNITS OF KSI. 

X5=X4/Y1 

C X5 IS THE PtASTICTTY CONDITION PATIO BASEO ON MASUBUCHHS STRESS 

C INVARIANT. WHEN IT IS E.T. OR G.T. i.O PLASTIC STRAINING FAS 

0 OCCURRED. 


RETURN 

END 
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WELOtNO STRAIN ROSETTE CAT A! EXPER I RENTAL RESULTS 

• 9 * 6 M*#**********#*********. »*■«••«»»••••••••*••«« ••**»••*•»■• 


JEST NUMBER Si BUTT VELD 

ALUMINUM UlPT 60 b 1 IN T*41 CGNOITIPN 
THICKNESS « 0*250 INCHES 
CM HELD PROCESS 
ARC VOLTAGE » IN. 

AMPERES - 2S0. 

SPEED OF TRAVEL • 32.16 INCHES /HI NOTE 
6063 FILLER METAL* 0.0625 INCH CIAHETER 
UNITS OF STRAIN * HICRCSTRATN 


— samsi. 

TERR ■ 0.00 rCRPuTATICNS REPRESENT 
BLANK SPACES ON INPUT CAROS and are 
MEANINGLESS. IGNORE ALL SUCH COMPUTATIONS. 
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CAUTILN 
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UANSTtRSt tillTANtl ♦ - “ »M J IN U‘i> 

TIN* I.OU 1.14 

( > EK’N-US ) CwlfJTt«iINtI (n.J I * /• l I L- «■ 1 


0.00 


TtMPERATuAt tut (..PAH* 1 l». J3 

7 6. j. 

STRAIN A. IMlaSI 

U.JU 

Vt.xj 

STRAIN H. (ALAS I 

J. 1 J 

IV. JJ 

STRAIN c. (Ml -SI 

0 4. JJ 

9.UU 

STRAIN a* IHlLrt* 

1.8V 

bt>. i r 

STRAIN *J. Mlt«» 

t.sv 

20. j* 

STRAIN C. (MtLMi 

»6.dV 

lliil> 

STRAIN X. WUrtl 

•6. *9 

s J. y» 

STRAIN T. I*LIHJ 

1.6V 

20. 1- 

GAMMA XT. ( Ml CHI 

-6J. JO 

50,4. 

PNIN.stR.l. (- 1 

7 7. >v 

t-S.O' 

PKlN.STR.2. 1 • 1 

-U.lt, 

<-./ i 

PHI A (i}LURfct4* 

XX. >J 

29.42 

SI ,MA XA U>1 ) 

ii.r* 

G.bS 

SluMA TT U 4 II 

0.2 7 

tl.6 > 

TA’i at Ins I i 

-0.2* 

0.16 

TItLJ SThL44 IKSII 

ii.r* 

IV.’IV 

t USD 

J.S4«9t IJ 

J.IASjt JJ 

I/AIEU S1ROS 

O.UmBE-OI 

Q .2 >86 1 - 01 



' _ , u.jj 


TCNPEXATUREIDFG.MHPI 

70.110 

78.00 

*4.0-. 

TE APEMTjktlutu.E-H*! 

r».uo 

7b. UO 

STRAIN A* |MEA$) 

-10.00 

-84.01 

STRAIN A. (NtA,l 

-18.00 

1 l.uO 

srittN a. iheasj 

32.00 

-39.00 

STRAIN 8. (MfcAS! 

lu.'JJ 

14. VO 

STRAIN C. (NEASI 

54.00 

-1,00 

STRAIN C. INEAsJ 

-/7.jo 

-13*00 

STRAIN A, (NECMI 

1.89 

-to . 11 

STRAIN A. IMbl.*** 

-l».(l 

-46. II 

STRAIN 0. (NECHI 

37.8? 

73.89 

STRAIN 8. (Nbi.NI 

Sl .84 

1.89 

strain C. (NECK* 

15.39 

-10.11 

STRAIN C. (HU** 

-84.11 

-22.(1 

STRAIN X. (NECmI 

-16.11 

-94,11 

STRAIN X. (NECHI 

-iso. 11 

•TO. 11 

STRAIN T> (NECMI 

’7.89 

73.89 

STRAIN Y. (Htcnl 

sl.o v 

1.64 

GAMA XV. (NECHI 

-18. JO 

C.J3 

GAMMA AT. IMECHl 

7/. JJ 

-24.33 

PafN.STR.l. I s I 

39.15 

73.85 

PklN.STK.l. I " 1 

A9.ZB 

3.84 

PRIN.STR.Z. t* 1 

-17.57 

-94.11 

PKIN.STm.2. (* 1 

-143. SO 

-72.06 

PHI X (DECREES 1 

80.78 

90.00 

Prll X IDbOkblsI 

78.40 

80. 78 

SION* XX (RSI* 

-0.04 

-0.78 

SIGMA XX (AMI 

-1.42 

-U.7J 

SICNA YY (RSI 1 

C. 37 

0.48 

SIGMA VY OolJ 

- J.13 

-3.24 

TAU xy irsir 

-0.07 

o.oc 

TAU AY (Rill 

0.27 

-6.09 

VIElC STRESS (RSI* 

19.99 

3S.54 

YILLo SIXEII IRSII 

J9.49 

19.49 

1 IRSII 0* 

2292E DO 

0. UC5F 01 

I INS 1 1 

0.102Jt 01 

0.4607b OO 

I/VIELO STRESS 0. 

5712E-02 

0.2763E-01 

1 /YIELD SlKtss 

J.4oSl£-Jl 

0.1 1 37 1- Jl 

15.00 

tfnp er a tur e i de o. r aha i 

78,00 


15.00 

TENPEPAT0RE(Ofct».E ArtRl 

78.011 

79.00 

strain a. ineasi 

-10.00 


STRAIN a. (NEASI 

6. JO 

-7.0.0 

strain a. (NEAS* 

18.00 


STRAIN 8. INCAS! 

54. JJ 

11. JJ 

STRAIN C. INEASI 

-84. 00 


STRAIN C. INEASI 

-123.00 

—15.00 

STRAIN A. (NECMI 

1.89 


STRAIN A. (HUH! 

7.89 

-64. J2 

strain a. (NECMI 

43.89 


STRAIN 0. IHtUlf 

35.89 

19.98 

STRAIN C. (NECHI 

-118.11 


STRAIN C. (Nttrtl 

-184.11 

-22.02 

STRAIN. X. (NECHI 

- 160.1 1 


STRAIN X. 1 MICH! 

-232.11 

-1J6.J2 

STRAIN V. (NFCHI 

43.89 


STRAIN Y. INUK! 

55.69 

19.4b 

GAMMA XV. (NECHI 

120.00 


GAMMA XT. (MtLMI 

192.00 

-42. JJ 

RRlNfSIA.I. 1* 1 

60.23 


PRIN.STR.l. (•• 1 

•4.96 

IJ.JV 

ERIN.STR,?. (■ I 

-176.45 

. 

PRIN.SlR.2. (• 1 

-261.18 

-109.43 

PMI X (DECREES 1 

74. 7T 


PHI X ( DEGREES I 

73.13 

83.78 

SIGMA XX (KSII 

-1*64 


SIGMA XX IRSII 

-2.41 

-0.43 

SIGMA VY IRS I) 

-0.11 


$ 1 uHA YY IRSII 

-0. 24 

3. J2 

TAU XY IRSII 

C* 45 


TAU XY IRSII 

0.72 

-0.16 

YtElO STRESS (ASTI 

39.99 


YIELD STRESS IKSIJ 

39.14 

*9.99 

i (RStl 0. 

1971E 01 


I IRSII 

J. 2128b J1 

0.6365b 00 

I/TIELO STRESS 0. 

4928E-JI 


I /YIELD STRESS 

0.682 It-Jl 

0.1592E-J1 

20.00 

TE“PFRATOR.f lOEG.EAHR I 

78.00 

78.00 

14.00 

TENkfcPATURtllUV.EANR 1 

■ to . JO 

79. uO 

STRAIN A. (NEASI 

38.03 

94.00 

STRAIN A. INEASI 

24.00 

[07.00 

STRAIN R. (NEASI 

90.00 

43.00 

STRAIN B. (NLASl 

6 J. JJ 

3 7. 33 

STRAIN C. (NEASI 

-204.00 

-43.00 

STRAIN C. (MLA*i 

-Jls. JO 

-87.00 

STRAIN A. (KECH» 

49.89 

169.85 

STRAIN A. IMtCMI 

25.-9 

49.48 

STRAIN B. (NECHI 

53.89 

175.84 

STRAIN M. (MfcCHl 

61.84 

19.98 

strain e. (nechi 

-238.11 

-52.11 

SIkAIN C. ffU.HI 

-376,11 

-96.02 

STRAIN x* (NECHI 

-244.11 

-58.11 

STRAIN X. (MtCHl 

-412.11 

-64. J2 


152 



STRAIN Y. t"FCM) 

ss.it 

T75.B4 

c*w», XT. tKECHI 

IBB. 00 

222.00 

P«lN.ST«.l. t* ) 

I 13. B? 

220.IT 

PPIN.STA.2. 1- 1 

-102.04 

-102.39 

PHI X IDEGRfFSI 

66, OB 

6B.25 

SIGH* XT USD 

-2.54 

0.01 

sigh* rr txsi> 

-0.29 

1.77 

HU XT USD 

1.08 

0.94 

YIELD SMfSS IASI) 

19.99 

39.99 

r cash 

0.3012E 01 

0.2369E 01 

1/VlFlO STRESS 

3.T511E-31 

J.5924E-Q 1 

24,<3j 



Tf MPENATUR El OEG.FAN** TS.OO 


STRAIN A. CMEASI 

ITO.OO 


STRAIN 8. IMFAS1 

14.00 


STP AIM C. IMEASI 

-454. 00 


STRAIN A. TMFCH) 

101.89 


STRAIN B. INECHI 

19.09 


STRAIN C. I*1£CMJ 

-490.11 


STRAIN X. IMECHT 

-128* 1 1 


STRAIN T. CNEO) 

19.49 


GAMA XT. TMECHJ 

672.31 


PRIN.STP.t* I" 1 

224.27 


PRIN.5TP.2. I- ) 

-512.49 


PHI X (DEGREES) 

5B.69 


SIGMA XX fXSSI 

-3.61 


SIGMA YY |ASI» 

-l.Ol 


TAU XT (RSI) 

2. 51 


YIELD STRESS IASI) 

19,99 


I IXStl 

0.4751F >1 


t/riEifl stress 

0.1064E 00 



25.00 


tfmperatupeideg.eamri 

90.30 

STRAIN A. I Mf AS J 

200.00 

STRAIN R. |ME as > 

56.00 

STRAIN C. INEAS1 

-408.00 

STRAIN A. IMFCMI 

212.09 

STRAIN B. IMFCMI 

62.09 

STRAIN C. TNECMT 

-441.91 

STRAIN X. WFCHj 

“191.91 

STRAIN V. IMECNI 

62.09 

GAKM* XT. IMECHI 

654.00 

PRIN.STP.l. «•* 1 

254.92 

PRIN.STR.2, ( • ) 

-404.74 

PHI X TOFGPEFSI 

59.21 

SIGMA XX (RSU 

-2.67 

SIGMA VV ISSN 

-0.0 1 

T Ai( XT IASI) 

2.44 

YIELD STRESS USM 

19.98 

I (RSI) 0. 

3BC6F 01 




strain Y. UUrtt 

el. 44 

19.94 

GAMMA AY. <Mti.H) 

YMi)J 

144. JO 

l*»I A. SIR. i. I" ) 

U/.uS 

61. J4 

PmTN.ST A. 4. 1" 1 

-S«Ji J t 

-I 35. 17 

PHI X (OEGnttSl 

69.49 

60.10 

SIGMA XX I4SI) 

—4.44 

-0.45 

sigma rr iasi I 

-0.55 

0. Id 

TAJ XT (ASH 

1.91 

0.54 

YIELD STmLaS (as) ) 

39. YN 

JV. 9V 

t IASI) 

01 

J. 13541 Jl 

1/rlcLJ sIhls* 

0.1. 46c. JO 

v. 16 A6f - Jl 


41.30 

TtNPE RATcjAe lot G.EAH, J 

lA.tfO 

74.00 

SIXAIN A. <Ml-j) 

74.00 

4 7.00 

strain 6. (ME A. t 

60 .JO 

25.00 

STRAIN C. INCAS) 

-347. Jo 

- J09.o0 

STRAIN A. (Mitrt) 

7 3.49 

-l J. 34 

STRAIN fl. 1 HU. HI 

ol.dV 

7.93 

STRAIN c. IMU.H) 

-jeS.U . 

- Jit. 04 

STRAIN X. 1MACHI 

-J7ci.il 

- : J -. a 

STRAIN Y. iMELrt) 

61.49 

7. No 

GAMMA XX. (MICH 1 

404. JJ 

JJG. JJ 

PRIN.SIR.1- 1- 1 

ibi .2u 

Lb*- — 

PKIN.STR. 4. (* ) 


— 3V —4 

PHI X (UtCREtol 

66. 74 

64. ;v 

SI-.MA AX IRSl) 

-4.01 

- J. J4 

SIGMA YY IASI 1 

-J. 71 

-3.47 

TAJ XY IASI) 

1.74 

).!' 

YIELD STRESS 1AS11 

J4.49 

14, 4 » 

I IASI) 0. 

4J>11 01 

0 . It 1st > 1 

1 /YIELD SIHESS 0* 

10431 oo 

0.91 **i( - >i 


2 3.00 


T £ toe HtURt CUty.l-AiM J 

XM.JJ 

74. J*. 

STRAIN A. (NcAI) 

A Oj- Ju 

1 '!•%) 2 

STRAIN B. 

JO.JJ 

Oi JJ 

STtt»IN C. IMt-o) 

-4oA«0u 

-44 (. ij 

STRAIN A. Intt-Hl 

169. Jo 

1 1 ' 

SMH 4. (MtCHj 

Jl.94 


STRAIN c. I Ml .ft) 

-SI'** Ja. 

6 j/ 

STRAIN X. MtCn) 

- 3 76, J4. 

“if . it 

strain r. (MitHi 

Ji mfO 

*• 

otMRA XX. Mi CHI 

tl OS » 4/ v 

V. ** «■ 

P4IN. STR.i. | 

6. aIC* * d or 

l O M t / * 

PRIM. SIR. 4. I" 1 

-a 7 J.4- 

-4'-4. 

PHI X IL'tGtct.I 

U J* H l 

4 7. j f 

S!o«A ** IXilJ 


* >.*‘J 

Sl.RA t Y U»U 

- 3. Jt, 

-1. i 7 

III XV (Rill 

2. >7 

i ' # 4 1 

TILL J ST At So CN.ll 

39. 

r-.Vr 

1 (Xill J 

l Jl 

J* 1t*lF ji 


153 


OBWMABM 
OP POOR flP 



l/VI ELD STRESS 

0.9319E-01 

26.00 

TeNRFRATUAf (DFG.MMR 1 

i 93.00 

STRAIN A. (REAS) 

302.00 

STRAIN A. INCAS) 

3*8.00 

STRAIN C. INEAS) 

•266.03 

STRAIN A. INECKI 

316.56 

STRAIN P. INECHI 

3T6.S6 

STRAIN C. CNECM) 

•270.66 

STRAIN X. INECMI 

-339.66 

STRAIN V. INECHI 

376.56 

GANN* AY. CHECH I 

596.00 

RRIN.STR.t. !• ) 

691.93 

PRIN.StR.2. »• 1 

•666.99 

RM1 x (OECREES) 

70.12 

SIGN* xx IASI) 

-1.65 

SIGN* TV IRSII 

3*9t 

TAU XV (RSI) 

2.23 

YIELD STRESS. IXSII 

39.9* 

1 (RSI) 

0.3658E 01 

1 /YIELD STRESS 

0.1366E 00 


26.5 J 


T f NPF« ATu» 1 1 CE G. RANR 1 

98.00 

STRAIN A. INEAS) 

6*9.00 

STRAIN A. INtASI 

986.00 

STRAIN C. INEAS) 

•266.00 

ST'-AtN A. INECN) 

5)6.93 

STRAIN 9. INECHI 

596.93 

STRAIN C. INECHI 

-293.07 

STRAIN X. INECHI 

-383. )7 

STRAIN V. INECHI 

596.91 

GANNA AY, INECHI 

798.00 

RRIN.STR.t. 1- 1 

737.06 

PP|**.5tR.2. |- | 

•323.20 

PHI * infGRffSI 

70.39 

SIGNA XX IXSt) 

1.79 

St&NA YY |K$I 1 

9.11 

Tit! XT IRS II 

2.99 

VtEtO STRfSS CASH 

39.86 

T US!) 3. 

9690E 01 

I/TIELC STRESS 0. 

222IE 00 



27.00 

TENPfRATUREIOFG.EiMa J 103.00 
STRAIN A. <*fASt 630.00 

STRAIN *. INfASI NS8.0D 

STRAIN C. INCAS) 

STRAIN A. |NfCHl 6*9.0* 


I/YtElO ilKLii 


».1UA JO 


3.96 07fc- Ji 



26.)) 


. 


TtNPcRATUKEI 01 G.FArtrf 1 

79.00 

79.00 

STRAIN A. 

INtASI 

21). 3) 

213.)) 

STR At K ». 

INtASI 

60.00 

37.00 

STkAIN c. 

INEAS) 

•6)7.00 

-605.00 

STRAIN A. 

INtlHI 

211. OR 

157.96 

SIKAIN 9. 

INtlMI 

61.98 

19.98 

STRAIN C. 

INtLM) 

— 6 IR. J2 

- . , 2. 02 

SIKAIN X. 

(NELH) 

-3tV.32 

-276.02 

STRAIN V. 

INtLHI 

61.96 

19.98 

GANNA XT. 

INtCHJ 

690.00 

570,00 

PA1 it. SIR. 

1. I" 1 

263.26 

19). 66 

PmIN.STR.2. I" 1 

-329. Jl 

-66 7.69 

PH] a IDEGRitS) 

39.16 

56.o6 

SIGNA XX 

(K3D 

-3.27 

— 2.r2 

SIGNA Vt 

IXSi J 

—0 . )6 

-).6l 

Uk> XV 

IXSII 

2.39 

2.16 

YtEL.) STRESS IKStl 

*9.99 

39.99 

1 (RSI) 


O.Ntoit Jl 

V. 36111 Jl 

1/fttLJ 3lkt»3 

9.10661 00 

0.90311- Jl 


26.5) 


T E 0 Pt RAluK t (UtG.T AHA ) 

79.0) 

80.0) 

STRAIN A. 

I ib*S) 

2nd. 00 

<65.00 

STRAIN 3. 

INLA3) 

166. 00 

165,00 

STRAIN C. 

iRt*3) 

-3)9.)) 

-605.0O 

SlRAlN 3. 

1 iUn» 

289.98 

188.09 

STnAIN u . 


163. 98 

126. 09 

Sit*;* c. 

INLLHl 5 

—600 .02 

-611.51 

strain x. 

I Xt.HI 

-230. 02 

-351.91 

SI tAIN T. 

iNttm 

163. ■*» 

126.09 

GiHNA XT. 

IRiLHl 

uVU.OJ 

JO 

RkIN.STM. 

1. I* ) 

366.26 

2)2.28 

PdlS.StR.*, 1* 1 

-636. jil 

-696.99 

PHI X tOtuKlt*) 

60. 1 L 

66. *) 

SIGNA XX 

1X31) 

-2.13 

- 3.1U 

SIj"V TT 

(Ait 1 

0.88 

6.31 

TAJ XT 

tX3l ) 

2.89 

2.26 

Tit Lo i 1x133 (R31) 

39.99 

39.96 

I (KS1I 


0.3d /RE Jl 

3.6268* 01 

I/TItLO 311*633 

0.96991- )1 

J.ITtTt Ou 


23.00 

T c 4(*f RATuKt lOtU.R AN< I 63.03 9J.J0 

ST. AIK A. IJU.A3) 286.J0 

STRAIN rt. INtASI JId.Ju 

STRAIN C. INtxSl J) 

STRAIN A. tNUH) *90. 56 It !.•><* 


15 A 



STRAIN 9 . IXECH) 

651.38 

STRAIN C. INECHI 

- 410.91 

STRAIN X. (HECHl 

- 392.91 

STRAIN T. (NECHI 

651.08 

GANNA XT* (NECHI 

toso.oo 

MtN.SW.t. 1 " 1 

880.14 

PRIN.SW. 7 . I* ) 

- 621 . 9 ? 

•HI X (DEGREES) 

6 T. 0 I 

SIGMA XX |k$ t 1 

2 . 86 

SIGMA TT IKS II 

10.66 

TAU KT IKS! I 

4.13 

Tlfto STRESS (KSU 

34,82 

1 WSI) 

O.IOITE 02 

I/TIELD STRESS 

0 . 2 SSSE 00 


2 7.50 


TERRE* ATuREtOIC.EAHRI 

I 2 S.OO 

SWAIN A. IHEASI 

570.00 

SWAIN 1 . WEAS) 

444,00 

STRAIN {. IMFASI 

- 432 . JJ 

STRAIN A. IMECH) 

609.49 

STRAIN 9 . INECHI 

519.49 

STRAIN C. INECHI 

- 416 .S 1 

STRAW X. (NfCMl 

-S 26.51 

STRAIN V. (NECHI 

514.49 

GAnma XT. (NECM) 

1026.00 

P 8 IN.SW. 1 . 1 " I 

761.39 

r: IN.STR. 2 . (• I 

- 560.41 

RH| X I DEGREES 1 

64.15 

SIGMA XX IKSII 

T .39 

SIGMA TV IKSII 

13.66 

T*U XT IKSII 

3.91 

TlEtO STRESS (KSU 

39.64 

I (KSTi 

0 . 1232 E 02 

UrlElO STRESS 

0 . 310 TE 00 


29.00 


T EMREA ATUR § | oe G. f AMR I 

154.00 

STRAIN A. IMFASI . 

194.00 

STRAIN A. (MfASI 

350 . JC 

STRAIN C. IHEASI 

- 264.00 

STRAIN A, (RECHI 

245.50 

STRAIN 0 . IMECH) 

345.50 

STRAIN C. (MfCHI 

- 258.50 

STRAIN X. IMEChI 

- 400.53 

STRAIN T. (NECHI 

395.50 

GANNA XT. WECH) 

504.00 

RR 1 N.SW.E. !• | 

467.95 

PRIN.SW. 2 * (• 1 

- 400.46 

•HI X IOCGREESI 

73.96 

SIGMA XX IKSII 

11.62 

SIGMA TV USD 

IT. 5 ) 


SMIU 9- INttnl 
SWAIN C. («CCHJ 
SWAIN X. (MttHf 
STRAIN V. I MU. HI 
GANN* XT. IMtLH) 
PKIN.STR.l. I" I 
PH 1 N.SW. 2 . t" 1 
PHI < lOLMUtl 
Slut* xx 1**0 
Slit* TV IK *11 

r*j xt tx>u 
VltlJ >Wfca* Wall 
I (KSU J 

1/rIEtJ iWfci* 0 


25 . 5 u 

T EWt:»*TuRclLitv.t-A>M3 
SW 4 IK t. (Kt-lJ 
SWAIN 4 . 

swaik c. (iu-si 
SWAlN A. Wt-HI 
SWAIA A. (Nt-«l 
SWAIN C. Ml-lll 
SW.il I X. 1 ML-H) 
SWAIN T. (MtLHJ 
GARM* XT. IHtu'tJ 
P- 1 N.SW. 1 , |“ ) 
PHh.iis.i. f" I 
PHI t (OcunetSl 

sjut* xx ikxii 
S lut* TT (KiU 
TAJ XT IKSIJ 
VIlIJ MkEoo Will 

i w>n 

l/tltllJ OWtSi 


26.00 

T t 4 PE«A 1 u*tluCip.E*m) 
SWAIN A. INUTI 
SWAIN 8 . IXt.il 
STRAIN t. (iUA*J 
SWAIN A. I tti.n| 
SWAIN 8 . (RtLHI 
SWAIN C. IMtLnl 
SWAIN X. tNti.nl 
SWAIN T. INCUR) 
(•ANNA AT. IMU.nl 
P-IM.SW.l. (• ) 
•MV.SW. 2 . I” t 
PHI X IDtUtttil 
SIGMA XX IKS 1 1 
SIGMA TT (KSU 


4X3.64 

. 376.54 

-409.46 

-255.16 

-339.46 

-339.-6 

320.64 

326.54 

600.00 

524. uj 

446.52 

421.44 

-4*5.45 

-404.77 

60*06 

n. is 

-1.65 

*1.29 

4.40 

3. —o 

2.25 

1.94 

49.NO 

*9. MS 

5 JbOt 31 

0.44271 J1 

12 lit 00 

0.17343 -jJ 


oV. Jj 

64.33 

4 54.00 

52 4.00 

522.00 

54 4. j _ 

■-X44.0G 

- 719 . 01 * 

4l t. JO 

44 7.6 6 

525.46 

56 7.00 

— 302.12 

-244.12 

-406.12 

-434. 12 

525.40 

567. 40 

7 64.00 

7 32. uJ 

654. 47 

69 3.26 

-510.61 

-526.57 

70.-0 

71.51 

-0.2- 

-0. 2V 

- 6.60 

7 ■ ju 

2.79 

2. 7 J 

49.92 

49.N2 

74101 ol 

4. 7 rest 41 

104 It Jj 

0. In 331 no 


1U4.UO 134. JJ 

046.00 TOT.oJ 

t> 06 . JG 649.00 

-351.30 -417. JO 

»nS,A4 651.55 

615.04 ' €JN. 55 

>>1 -4W.15 

•ir-.ii *i , .d..s 

ois.ju e>v. 55 

i-SO.OJ 1074. jj 

44 1.65 ttr. Jt 

Hwl.Tt ~t.2c.it, 

H>.AA til.ul 

2 .N 3 2 .n 5 

10.43 io. n 


155 



T4U *r 

mn 

I. as 

YIELD STRESS f«Slt 

)«.!« 

i i«$t» 


J. 15b If 32 

t/VIFLD 

sm» 

0.39736 00 


29.00 


TEnpERATUBECDEG.FAhA I 

219.00 

STRAIN A. MFASl 

-526.00 

STRAIN ft. IKMSI 

470.00 

STRAIN C. t*F4SJ 

-66.00 

STRAIN 4. IHffHl 

-441.64 

STRAIN Q. INfctMl 

548.36 

STRAIN ,C. IHtCHI 

-27.64 

STRAIN X. lUfCrtl 

-1017.64 

STRAIN V. 1 4 EC HI 

648.36 

GAMMA IV. IHfCHl 

-414.00 

P* IN .STB. | . 1- 1 

575.26 

PPtN.SfB.2. 1* 1 

-1044.54 

»Ht K lOfGRfFSl 

62.60 

SIGMA XX <K$I I 

17.24 

SIGMA VV (HSU 

28.48 

Tfj xt ms ii 

-1.49 

YIFLO STRESS USD 

38.40 

i <*sn 

0.7476E 02 

mvtfto STRESS 

0.64476 }3 


30.ro 

TFMPEAATUBEIOEC.fAHBI 

260.00 

>5.00 

strain a. <«as» 

-866.03 

240.00 

STRAIN ft. IHfASI 

674.00 

237. OC 

STRAIN C. <*EAS» 

46.00 

-61.00 

STRAIN ** INFCHl 

-785.64 

313.89 

STRAIN ft. (HECH! 

750.36 

349.69 

STRAIN C. INKHl 

64.36 

-70.11 

STRAIN 1. (*ECH| 

-‘451. 64 

-106.11 

STRAIN V. iNfCHJ 

750*3, ■ 

349.89 

GAMMA *V. 1*ECHJ 

—BTC. 33 

384.33 

PRlN.STR.lt f* 1 

833.18 

419,96 

PAl*.ST«.?. (• 1 

-1534.46 

-176.16 

PHI X (Ofc&Rf f S 1 

79.22 

69.95 

SIGMA XX 1*511 

21. HI 

0.12 

SIGH* TV IRSII 

37.28 

3.55 

TAO XY 3 BSD 

-3.06 

1.44 

YIELD STRESS IXSII 

37.32 

39.59 

1 tKStl 

0.3230E 02 

0.4064F 01 

I/YIELO STRESS 

0.6655E 00 

0.1016E OO 


11.00 

TFNPERATuRFICtG.fAHRl 28*. J3 
STRAIN A. IMMSI -1018.00 


760 XV 

IKSU 

3.92 

4.01 

YIELD STRESS IRill 

39.62 

3V.dt 

1 (Kill 


3.963 7t Hi 

3.1C19E 32 

1/VIEL0 

SlHfcSS 

o.247oE 00 

0.2559E 00 


2t.50 




T tHPE SATURfc IGiG, b AH* 1 

129.33 

129.00 

sum r. 

IHtASl 

654. JO 

713.00 

STRAIN 8. 

1 NLAS1 

436. 33 

487. 33 

STmAIN c. 

ixmi 

-441.30 

-483.30 

Sr-tAI ft A. 

tNEl-rtJ 

6 7 4. 1 » 

678.11 

STRAIN ft. 

IMtCnl 

403.11 

492.11 

STRAIN C. 

tHECHl 

—4 79.0V 

-467.09 

STRAIN X. 

(HtCuJ 

-201.09 

-281.04 

SIXAIN T. 

i.NtCHt 

400*11 

44.. n 

GANNA XY. 

(NLLH) 

1150. JO 

1 146.00 

P3IN.STB. 

1. 1“ 1 

742.22 

716.56 

PB 1 N. SIM.X. 1* 1 

-344. 3 J 

-50b. 33 

PHI X lOtbNLtjl 

ol. t>7 

62.02 

StoNl XX 

IMI 

6.51 

0.55 

slitHA rv 

IKSU 

14. Lb 

14.29 

TAJ XV 

(Ail 1 

4.24 

4.25 

VtcLi) STkES> INSlI 

39. uO 

34. t.3 

1 IKSII 


3.1265E 32 

3.1x46b 32 

I/TIELO sIhESS 

3.3240b 00 

0.3. TIE 30 


2 7.00 


TEIPEPRToKEIUbu.bAHAl 

lo5» 3 J 

16 3. 33 

STRAIN A. IHtASl 

342.03 

341.30 

SIXAIN 0. 4Hb Ajl 

234. 4|0 

307.00 

STRAIN C. IHtASl 

-334.00 

- l63»UO 

SlHAlN 4. lMU.Nl 

391.09 

329.73 

STRAIN 8. ( AfcCri 1 

203.39 

335.73 

STRAIN C. t MtLHl 

-352.41 

-324.27 

STRAIN X, INtLrtl 

-244.41 

-33G.2I 

STRAIN V. (NbLnl 

263.09 

335.71 

GAMMA XV. I HUH I 

744. m) - 

654.00 

PhIN.sIB.I. (■ 1 

415.26 

46 9.44 

PAIN. SIR. 2. 1- ) 

-437.37 

-46 1.48 

phi x loe&Actsi 

62*60 

b 7, 7t 

SIGMA AX IXiil 

15.13 

I 9.00 

S11H1 YV IKSII 

18. it 

u. «a 

tiu XY (Kill 

2.72 

2.43 

YIELD sTkLSS IKSII 

34.21 

39. ti 

1 [Kill 

0.1762t 02 

D.1719F 32 

1 /YIELD STRESS 

3.4494C 00 

0.43607 03 


27.50 



T LNPtBATOktt Oto.bAHX I 

201.00 

198. GO 

STRAIN A. 1 HE AS 7 

-96.30 

-41.00 


ST* At* ft. IREKI 
* STRAIN C. INMSl 
STRAIN A. IHECHI 
STRAIN e. IHECHI 
STRAIN C. IHECHI 
STRAIN A* IHECHI 
STRAIN V. WitCHI 
GANNA AT. IHECHI 
PRlN.STR.l. r I 
PRIN.STR.2. (■ I 
PHJ X I DEGREE SI 

signa xx txsn 
sicna vv iksi i 

TAU XV (KS1I 
VIELO STRESS IXS1I 

l ixsn 

I /VI ELD STRESS 


T7B.00 
-24.00 
-931. 5fc 

938.34 
-3.66 

-1793.86 

838.34 
-948.32 

971.10 
-1876.4? 
80.13 
23.35 
41.59 
-3.28 
36.5Q 
1.JMIE 02 
0.9854E 00 



33.00 

TENPEMTURfl.OEG.EAMRI 305.30 
STRAIN A. IHEASI -1042.00 
STRAIN 8. (NE AS) 848.00 
STRAIN C. INEASt -144.00 
STRAIN A. IHECHI -995.09 
STRAIN 8. INECM) *98.91 
STRAIN C. IHECHI —l* 3.09 
STRAIN X. IHECM1 -202?.?9 
STRAIN V* IHECHI 888.91 
GANNA XV. IHECHI -852.00 
PRlN.StR.I* (• I 949.8? 
r'RlN.STR.2. 1“ I -2088.03 
PHI X (DEGREES ) <1.96 
SIGNA IX IKSI I 24.58 
SIGNA TV IKSI I 44.57 
TAJ Xt IK5II “2.92 
VIELO STRESS IKSII 33.81 
t IKSI) 0.3R55E 02 
I /VI ECO STRESS 0.13TTE 01 


33.00 


1 € H£> ER ATUR E ( DEG. P AHRI 

309.00 

STRAIN A. 

IKEAS) 

-9T6. 19 

STRAIN *. 

IHEASI 

818. 00 

STRAIN C. 

IHEASI 

-284.00 

STRAIN A. 

IHECHI 

-934.04 

STRAIN 8. 

IHECHI 

853.96 

STRAIN C. 

IHECHI 

-288.14 

STRAIN X. 

IHECHI 

-2076.04 

STRAIN V. 

IHECHI 

853.96 

GAKHA XV. 

IHECHI 

-646. )J 

PRlN.STR.l 

. I- 1 

869.15 

PRIN.STR.2 

. !• 1 

-2111.22 


STRAIN 9. IHtAil 

14«..3> 

* 323. 3J 

Gin. 41 K C. (HEaW 

-i.4 7.JU 

-201.03 

STRAIN A. IHtLrtl 


-81.60 

STRAIN ft. (MlLnl 

2bl. >1 

3/4.20 

SIXAIN C. (HtCH) 

-14J.4V 

-141.60 

STRAIN X. (Hi ( HI 

-429.44 

-59 7. 63 

STRAIN V. * INc.Hl 

261.31 

374.20 

GANNA XV. CMitni 

114. 33 

6 3.03 

PRlN.STR.l. 1* I 

266.1* 

375.12 

PHI H.STK.2. 1 " 1 

-413.1/ 

-596.73 

PHI X IDEGHttsI 

63.31 

88.23 

SIGNA XX (KS1I 

IV. 91 

l/.vt 

3 IGHA TV (Kill 

26.99 

2 '..44 

TAJ XV IKSI I 

3.41 

J.2X 

vtetJ siKts> (Niii 

39.66 

19.72 

1 (HSU 

3.2nd 4 t 32 

(J.2224E J? 

I /VIELO SlKESS 

0*>407l 30 

0.9/5/ t Jj 


28.00 

1 EHPtR ATUot lu tU.E Amt, 1 

235. Jj 

2/5.30 

STRAIN a. (MtASI 

-426. JU 

-373.03 

STRAIN 9. IHlAkl 

289. J J 

451. 1) 

STRAIN C. IMfcASI 

-33.30 

-99.00 

S TkAI N A. (NLLHl 

-346.08 

-155.54 

STRAIN 8. IHtCHI 

365.12 

5jN.n1 

STRAIN C. INtCnJ 

-16.69 

-31.54 

STRAIN X.. IHtCHI 

-732.69 

— fe45a 54 

STRAIN V. IHtCHI 

365*1. 

5Ja.nl 

GANNA XV. 1 NLCHI 

-330.30 

— 324. JO 

PRlN.STR.l. 1“ I 

3*9.36 

5/6. 66 

PRIV.SIR.n. «" 1 

-757.16 

-414.36 

PHI A tOEGNttSI 

81.66 

8 5.33 

SlUHA AX (XSU 

23.41 

20.33 

SIGNA VV IKSII 

31.22 

10. J) 

TAJ XV InSII 

-1.17 

- 1.16 

VI EL J STIES) INSIl 

3/.9f 

JU.14 

1 IKS 1 1 3. 

<,6361 J2 

3.2b? >1 j/ 

I/VIELO SllttSS 0. 

7394 fc 03 

0.6Wt 30 


23.30 

I tHPEKAl JKL l Jtv.E AHA) 
SliiMS A. IMfcAS) 
STRAIN B - (6C»S) 
STRAIN C. (HtA>l 
S lit At N A. Iltu.nl 
Sin AIN 8. INtCrtl 
SI MAIN C. IHlCnJ 
STRAIN X. tMtV.Nl 
STii At A V. IMtCHi 
GANNA XV. (Ntv.nl 
PN IN. SIR . 1. I» I 
PRIN. SIR. 2. I* I 


263.33 

246. Jj 

-63' 1 

-57 1.33 

63 

55 4.30 

3.JJ 

-M.jJ 

-564.66 

-55 J. 4/ 

> IV. 1 1 

tic. 43 

11.11 

4.53 

1062.94 

-1 165.4/ 

504.11 

tlt.5* 

—5 76. uo 

-55 a. 00 

563.21 

*39.14 

•1113.94 

-1206.13 


I- 

tn 


PHI X (DEGREES 1 

83.78 


PHI X IDKinllll 

IT.T4 

At. 11 

SIGH* XX IXSIt 

24.58 


SloMA XX I RSI 1 

2S.Nl 

2 1. J8 

stGHi ri msi» 

44.12 


SIGHA TV {KSI 1 

J6.64 

34. 30 

t*u *v ixsii 

-2.21 


I1J XY IKSO 

'2.02 

-1.4H 

YIELD STRESS IRSII 

35.65 


YItLJ SlRtSS IRSI 1 

3/.2i 

37. 79 

I (K$t) 

0.3871E 02 


1 IRStl 

0.J24 7t J2 

3.29671 J2 

I/VtELB STRESS 

0.10R6F 01 


1 /YtELl) SlxtaS 

0.67204 JO 

J.7bS9I JO 

43.0L 

TEHPEaATuREI DEDHAM* 1 

255.00 

78.00 

29.00 

TtNPEAATuRtlGto.ERtMl 2a o.uo 

262.00 

STRAIN A. IPEASI 

-838.00 

306.00 

STRAIN A. 1 Kfc Asl 

-774. uO 

-703.00 

STRAIN B. IMFASI 

61*. 00 

357.00 

S T a A I N 8. IMt ASI 

95<f. Jj 

6 78, 3 3 

STRAIN C. IHFAS1 

—384.00 

-81.00 

STRAIN C. INLASl 

4. JO 

-6 S. 50 

STRAIN A. (MECHJ 

-780.33 

378.89 

S In AIN A. IMtCHI 

■Ml.'i* 

-68J. .5 

STRAIN B. {NEC HI 

665.87 

468.89 

STRAIN rt. iHtAHl 

Jli.iT 

701.99 

STRAIN C. !hECn| 

-372.03 

-100.11 

STRAIN C. IMt-HI 

2. J* 

*(..45 

STRAIN X. INI CH| 

-ma. is 

-140.11 

STRAIN A. <Ht(,H| 

-1 324.66 

-1422.49 

strain y. imecmi 

'45.87 

468.88 

STRAIN Y, IHEcrtl 

bl4.j4 

713.55 

CSNMA XT. (HECH| 

-408.00 

680.00 

GAMMA XY. (HLGHl 

- 72 0. JO 

-64 .. y J 

PRIN.STR.l. I* I 

682.61 

547,93 

PRtN.STR.i. I" I 

6 I 0 • o6 

7ns. >9 

PRTN.StP.2. I« I 

-1834.67 

-268.15 

PAIN. ST*. A. 1* 1 

-liTf. i o 

- 1 4To. lb 

phi x I DEGREES 1 

85.34 

71.99 

PHt A lOEWUtSl 

79. J— 

<tl. 12 

SIGHA XX (KSI ) 

24.15 

-0.38 

SIGHA XX IrSII 

2 7.J4 

22.-2 

SIGH* Y¥ IRStl 

41.28 

4.58 

SIGHA YY IRSII 

4J. 8 7 

0 7. 9- 

T &li XV txsn 

-1.41 

1.81 

TAJ XY IRSII 

-2.49 

-7.-2 

YIELD STRESS IRStl 

36.18 

39.49 

Y 1 cLO STKtSS (HSU 

ili.JIt 

) 7.26 

I fXSIl 

0.35B6E 32 

3.5322E 91 

I 1 R> I 1 

0.354/E J2 

j. j2t i r 02 

T/YtELO STRESS 

0.991CF DO 

0. 1331 F 00 

^ /YI2U SlALsS 

0.48521 00 

J.6 7 51 F .0 

50.00 

TENPERATtIRfI DFG.*AHR| 

267.00 

77.37 

iO. JO 

rtHPEAAt JnEIOtG.PAtPl 

1 ju.jo 

286. * J 

STRAIN a. IHEAS) 

-778,00 

648.00 

STRAIN A. |NtA>l 

-o 7c. JO 

-12 i.Jj 

r’RAIN B. IHEAS 1 

410.00 

261.00 

STRAIN 8. 1 ml AS ) 

726.30 

HI.j) 

"AIN C. IHEAS! 

-396.30 

-493.34 

STRAIN C. IHtAs! 

-it J. JO 

-117. JJ 

STRAIN A. CHECH! 

-688.79 

721. P2 

SIkAIN A. ( HtiH 1 

-647.45 

- £25, j it 

STRAIN 8. IHECHI 

483*21 

373.82 

SIrAIN 6. IMcLrtl 

794. oS 

35*. 14 

STRAIN C. IHECHI 

-362.79 

-502.18 

STRAIN C. IMtLHI 

— 47.95 


STRAIN X. I«*ECH! 

-1544.14 

-154.18 

STRAIN X. IHLCll! 

-i 649.95 

- 1 749. 6 1, 

STRAIN V. IHECHI 

483.21 

373.82 

SIrAIN Y. IMLbHl 

794.o5 

094. J- 

GANNA XV, IHECHI 

-336.00 

1224.00 

GAMMA XV. IMLCMl 

- 790. uW 

- 756. J J 

PRIN.STR.I. i* 1 

497,03 

776.33 

PAlN.SIR.i. »* » 

olJ.Jo 

4 J 6*10 

PPIN.STR.2. I* 1 

-1558.61 

-556,70 

PR IN. SIR. 2. I* 1 

-175>. 4| 

-1oJ3.4. 

PHI X (DEGREES 1 

85.30 

56.67 

RMl X IOLGrcLS) 

81.90 

dl.’il 

SIGHA XX IRStl 

21.20 

-0.53 

SIGHA XX 1RS1* 

2b. 7i 

2 3. 9b 

SIGHA YV IRSTI 

35.40 

3.44 

SIGMA « IRStl 

45.44 

41. If 

TAU XV IRSTI 

-1.16 

4.61 

I All XT IRSII 

-2.9b 

-2.b2 

YIELD STRESS tKSIl 

37.12 

40.03 

VicL> slRtSs Kill 

35. 9 J 

36.50 

! (RSI 1 

0.3080F 02 

0.5271F 01 

I IRSII 

0.3V75E J2 

J. 36 30 7 02 

l /YIELD STRESS 

0.8298E 00 

0. 131 fif 00 

1 /YIELD STRESS 

0.U19L 01 

0.94441 


158 


TEMPEPATuAEIOEG.f AHR1 

77.00 

32. JJ 

TEMPERATURE (UEO.FAMt) 

328.00 

310.00 

STRAW A. IMEASI 

750.00 

STRAIN A. I MC AS 1 

-066. 30 

0.11 

STRAW 8. IMEASI 

255.00 

STRAIN 0. IMEASI 

«28.00 

0.00 

STRAIN C. IMEASI 

-528.00 

STAAIN C. IMEASI 

-225.00 

-217.00 

STRAIN A. IMFCHl 

825.82 

STRAIN A. IMECHI 

-661.66 

-28.32 

STRAIN A. IMECHI 

367.82 

STRAW 0. INtCHl 

032.36 

i 1.68 

STRAIN C. IMFCHl 

-538.11 

STRAIN C. iMtCHl 

-203.66 

-215,32 

STRAIN A. IMFCHl 

-82.11 

STRAIN X. IMECHI 

-1957.66 

-755.32 

STRAIN V. IMECHI 

367.82 

STRAIN V. IMICHI 

632. '6 

11.68 

GANNA XV. IMECHI 

1362,00 

GAMMA XT. IMECHI 

-550.33 

18 7.00 

PRIN.STP.l. !• 1 

860.02 

PR IN. SIR. i. 1" 1 

059*98 

61.16 

Pn|N.STR.2. t« t 

-576.39 

PR IN. SIR. 2. 1" 1 

-1985.27 

—286*01 

"Hi X IPEGREESI 

56.16 

PHI X (UEGREfcSl 

86.16 

72.50 

StGMA XX IASI I 

0.26 

SIGMA XX IKSU 

29.21 

60.63 

SIGMA VY 

3.65 

SIGMA VV IK .1 

68.07 

62. *7 

T AU XV CASH 

5.12 

T4U XV tKSfl 

-1.09 

0.06 

VlEtO STRESS IKSII 

60.03 

YIELD STREj. IK >1| 

36.06 

35.61 

I IKSII 

0.5272E 01 

1 IKSII 

0.6109E 02 

0.61611 02 

I / V I EL H STRESS 

o.mer oc 

1 /YIELD STRESS 

0.1202c ul 

0.1 169 E 01 

52.00 


35. JJ 



TFNRERATUREIOEG.FANR) 

78.00 

t tMPE PATURt lUtG.EAMKl 

321.U0 

311. 00 

STRAIN A. IMEASI 

882.00 

STRAIN A. IMLAbl 

-720. JO 

U.JO 

STRAIN A. IMEASI 

617.50 

STRAIN 8. IHtAjJ 

776. JO 

0. JO 

STRAIN C. (TEAS) 

-633.00 

STRAIN C. IMILASl 

-ITS. JO 

369.00 

STRAIN A. IMECHI 

855.88 

STRAIN A. IMEcHI 

-736.60 

-09.62 

STRAIN 8* IMECHI 

528.88 

' STRAIN 0. WEcKl 

709.32 

10.38 

STRAIN C. INECWl 

-662.11 

SIXAIN C. IMCtMl 

-622.0B 

- 16 8.62 

STRAIN X. IMECHI 

-16. U 

STRAIN X. IMECHI 

-1916.68 

- -00.62 

STRAIN V. IMECHI 

575.88 

STRAIN V* IMECHI 

709.32 

10.38 

GAMMA XV. I**ECH| 

1358.00 

GAMMA XT. IMLCHI 

-202. JO 

319.30 

PRIN.STR.I. !• l 

1007.31 

PKjN.SrR.l. {■ I 

79o.o6 

66. jO 

PRIN.STR.2. 1" 1 

-683.53 

PRtN. *jTR#2. 1" 1 

-1926. Ul 

-666.55 

PH I X IDEGRfESI 

55.67 

PHI < lUtCihCtSI 

07.03 

70.80 

SIGMA XX IRS 1 1 

1.11 

SIGMA XX IKSII 

20.2* 

38.65 

SIGMA VV IKS! 1 

5.82 

SIGMA VV IKSII 

66. 62 

62.17 

TAJ XV (RSI 1 

5.26 

Tad xy iksII 

-0.96 

1.09 

YIELD STRESS IKSII 

38.88 

YIELD SlKEs. I K SI 1 

J>»13 

33.37 

t IRSII 

0.6586F 01 

1 IRSH 

O. 6 J 06 E 02 

0.60821 J2 

I /VI El IT STRESS 

0.1666E 00 

I/TlELJ STRESS 

O.ll 36t Ul 

O.lISJt 01 

52.50 


60. JO 



TFMPERATURf IOEG.FAHRI 

80.00 

T tM°ER ATUht lufc L>. T AHR J 

1 £90. JO 

286.00 

STRAIN A. IMEASI 

1036.00 

SIR--. ' A. IMEASI 

-302.00 

0. JO 

STRAIN 5. IMrtSI 

615.00 

STRAIN 6. IMtAsI 

>66* JJ 

J. J J 

STRAIN C. IMEASI 

-385.23 

STRAIN C. IMEASi 

-6 29. JO 

-611. Jl) 

STRAIN A. |MECH| 

1112.09 

STRAIN A. IMtCHl 

-5 37.15 

-U..J3 

STRAIN 8. IMECHI 

728.09 

STRAIN 6. IMECHI 

606.05 

29.97 

STRAIN C. IMECHI 

-393.91 

STRAIN C. IMECHI 

-6*7. 15 

-371. U3 

STRAIN X. IMFCHl 

-8.51 

STRAIN A* IMLLM! 

-1593. IS 

-611. J3 

STFAIN V. ITCh) 

726.09 

STRAIN V. t Ht CHI 

oOti.oS 

29.9 7* 


159 



umi nr. 

1SOA.OO 

GARRA XT. (HtC.NI 

■ 

•90.00 

361.09 

W|K«SH.t. 1* 1 

1 ITT. AS 

PAIN. SIR. 1* (• 1 

609.7 A 

94.43 

MIN.UM. 1" 1 

-479*46 

PR1N.SI A.2. (■ 1 


-47S.4* 

Wt X lOUIEISI 

SC *05 

PHI X lOCGKUSl 

•4.6A 

70.35 

StGRA XX msn 

3*01 

SIGMA XX IKSII 

26. «2 

36.3S 

SICK* TV USD 

1.51 

SIGRA TV (KS1) 

41.96 

39.37 

TAU XV (XSII 

5*AA 

74U XT (Kill 

-0.il 

1.25 

vmo STRESS «X$I> 

».5( 

TIELO SlXtlS (XSII 

36. J 6 

36.22 

f CASH 

0.C5TAE 01 

1 (XSII 

0.363 IE J2 

0*3799 E 92 

f/VIEt.9 STRESS 

9.2I45E 09 

l/rttio STRESS 

0*102 OE 01 

0.1049E <11 

53.00 





TE KPER ATUR Et DEC. FAH9 I 

•4.09 




STRAIN A. (HfASf 

1254.09 




STRAIN 0. (REASI 

753*09 




STRAIN C. INEAT 1 

-409*00 

r 



STRAIN A. INK.) 

1321*72 




.STRAIN 8. (RECHI 

•60*72 




STRAIN C* (RECMI 

-417.24 




STRAIN X. (R(CH1 

44.72 




STRAIN V. (RECHI 

•66.72 




GARRA XT. (RECMI 

1746.09 




PRIN.STR*!* «• 1 

1420*63 




MtN.SfR.Z, !• 1 

-509.19 




PHI I (DECREES I 

57.61 




SI OR A XX IKSII 

4.91 




SIGH A TV (XSII 

11.09 




TAU XT (XSII 

6.56 




TIELO STRESS (XSII 

S9.9S 




t iKSII 

0.1060E 02 




I/TI ELD STRESS 

0.2652E 00 




53*50 





TEMP E R AtUR Et DFC.M MR I 

93.00 




STRAIN A. (REAS) 

1350.00 




STRAIN 1. (REASI 

70S.07 




STRAIN t. (REASI 

-39T.00 




STRAIN A* (RECMI 

1427.09 




STRAIN 5* (RECMI 

•21.09 




STRAIN C. (RECMI 

-402.91 




STRAIN X. (RECMI 

203.09 




STRAIN V. (RECMI 

•21.09 




GAMMA XT. (RCC«I 

1*30.00 




RRTN.STR.l. I* 1 

1477.66 




PRtN.ST*.2. (■ 1 

-453.67 




phi x (DECREES 1 

34,33 




SIGRA XX (ASM 

6.73 




SIGMA TT (RSI* 

12.89 




TAO XT (XSII 

6.06 




TIELC STRESS IKSII 

34.69 




I (RSI 1 

0.1218E 02 




I/VlflO STRESS 

1.39S4F 00 I 

■ 




160 


54,00 



TF*9E*4?U*ltDf<i.MH* t 


105.00 

4. ttffASI 


1 I TC.QO 

v»*n i. t«M$t 


403.00 

STRAIN C. t*fA$» 


-245.00 

UHJX *, 1 **£««> 


1252.03 

smai. •. i»khi 


404.03 

STRAIN C. U»KHI 


-245.4? 

IKECHt 


302.03 

ST*A1-, t, 


404,33 

Ktr- 


1519.00 

P«|N.ST9.1. 1* J 


1260.10 

•MS.SJK.J, |* i 


-2T4.04 

phi * mo*ersi 


45.16 

SKW XX IRStl 


11. T5 

<I$«a W l*S|l 


13.41 

TAU XT IhStl 


5.6? 

YtEtD STRESS InSI 1 


35. SO 

T IXS1I 


0.1332F 02 

t/rHUi stress 


3.33461 « 


5**55 ■ 


Tf VPEH4TU* M0FC.94H4 > ' 

123.00 

STRAIN i. WFASI 

rto.o? 

STRAIN >* (9IISI 

345.00 

STI»K C, !?'«» 

-145.00 

STRAIN A. (Fiu^l 

0T2.21 

STRAIN 9. (91(91 

4?6.2| 

strain r. (HccHi 

-135.T* 

ITIAh «« |MU9| 

240.21 

STRAIN 3. ICEtNl 

4T6.21 

&ANAA ??. IMECH1 ' : 

1308.03 

WI5.I1M. 1* 1 

833.65 

P«lM.SI«.2. 1* 1 

-I4T.23 

9H| A I0ECAFCS1 

51*05 

SIS* A IX IASI) 

13*3? 

Slolt n IASI? ■ ■ 

14.4? 

TlU *3 1*511 

3.T4 

YIHO STAFSS (A5t* 

36.64 

1 IASI 1 

J.1443E 02 

t/YlUO STAISS 

0*34589 00 


53.00 

Tf*ACAATUAE( 0 f 0 . 8 AMlt 

144*00 

STRAIN A. IDEAS I 

4 T 4.00 

STRAIN 6 . I«€A 5 I 

363.00 

STRAIN U 1 DFAS 1 

- 115*00 

STRAIN A. IDFCH 1 

5 * 0.45 

strain t. i*»tr*i 

5 ut .*5 




K 7 M 1 .V C. IttfM 



# -f*AI«# 1. ****** 

«*to*i 

-£»****. KT, f*FC*» 

P»f*»S<r». 1 . I* I 
I" I 

*«i * ictij^resi 
*!«**< mi » 
sigma tt ik*; 1 1 
taw *t iksii 
tirtc mu 

t usn 

jyri'i o s**f v* 


s*.oo 

ITMI» I. |l*f*H 

shaim n. imfasi 
witvr. i**»si 

STKAtM A. |*FCW 
SIMIN P. l«*K*l 
STP4t'« C. t*U*l 
». 1*100 
SI«AIN t. IRFfHl 
r,i»Nt n, lafCH* 
M>t»l.»f*.l. I» I 
wm.su.i. i- i 

wo * trtenrsi 

MW* K& IAS* I 
SIGMA TT mu 
TAW IT l«SII 
TIFtC S»«FSS mil 
I IKSII 

l/TI *t > ST*fSS 


6 T.JJ 

Tf 4PIK ATtfAS I*'! } 

ST*Ai,» i. 

suaim «* i*uu 
IIUIS C. I»*AM 

sf*Aiii a. (Mir nr , 

it«»H N tPECm 
Mim e. WffH) 
51 * 41 * 1 . IMltKl 
M»* 1 S V* «*«m» 
i v ,***«a : r*». i*Khi 
r if » 

msiST *.?. «« | 

■ phi x «t»te*H*i 

’• ." StfcMA *K WJtl : 
SICMA *X IKSII 
,TAU « IKSII 


... 



• II.M 
-16*55 
soj'.ak 

67?. 00 
671.76 

64. C< 
16.61 
IP-30 
2 . 4 1 
36 .*S 
1. 16621 7 ? 
0.42BSF 00 


166.0 ; 

276.00 
665.06 

5.IJ 

606.01 
6*7.01 

55.01 

-207.66 

66 T.it 

»5*.flC 

TCI . o * 

-?37.6> 

TP. 61 
20.61 
26.61 ' 
1.26 
18.6* 
0.266 IF 02 
9.62661 00 


213.07 

116.00 

563.00 
-t.OD 

256.83 

727.40 

S6.ee 

-612.23 

* 27.86 

204.00 
716.86 

- 621.25 

*6.63 , 
*4.0? 
32.21 
0 .M 




VKltt ”«t» tKSIt 
liitti 
l/viHit smss 


)l. 4 « 
O.MHI O? 
O.lIftF 00 


39.00 

iwusi 
i. «*c*9i 

i»*i» c« man 

«M W I* IMCHI 

v <**kh 

C. ««twl 
AT»*l** X. WICHf 
vm«* », i«fCHi 
«,»****.*¥. (MtHI . 
fH'tiJ'M* f * 
I**:.! 

Ml HOfCHfH 
iip**.»* t«$|i 

w ta$t I 

rw « »«$!» 

rlfct? H»f« KtU 
I l*5tl 

l/rlito iwm 


232.00 

* 2 .oo 

»«aj 

-7.00 
If 1.09 
7*3.03 

».oi 

7*3.05 

I32.0C 

W.« 

•9.95 
2 §.?S 
35.5* 
9.*7 
K .«4 
3.1210i 92 
O.MUE 09 


* 0.00 

UMHilTilMUftitlliOl 2 **.00 
1 *U« 1 . -IM .00 

$»Mt* k (DtMl 2 M .00 

Miraif t» fatal* : ,• -J 3 *.a 3 

*. f*€C**l ';. ;: L ‘ -«»*.*? 

siaaif o. mm )?«.» 

n»*m e. la^ct** * 2 M.*T 

traai* i. (fttcxt -V ■ ■ -lias.*? 

st*at* »* «*{•« , . ito.» 

<!«*** ■ tt. ■ -u*.m 

: 4.9 '.* . ItMl 

n I ».*?*. 2 . t* i -lut .*7 

Mi «: totofttsi JWl 

sieit* ii t**j* n.t» 

iiMl TV iiill ; : 1I*W 

litfiM- mil ■: *I.H 

*u*t $r*ts* mil IM» 

t mil ■ ■ IiIIME 92 

t/imo st»f« •♦fim m 


2*1.99 

-* 2.90 

S 79.00 

35.09 

tout* 

m.u 

93 . 9 * 
-3*9.0* 
>*!.** 
*.00 
1*1.9? 

- 9 **.l)« 
•«,«? 
31.39 
*0.70 
Ml 
37.29 
0.3*932 02 
0.99032 00 




271.90 

*300.99 

3*9.00 





227.09 

•ita.tr 

3*7.03 

-*24.91 
547.01 
-944.00 
195.31 
-973.27 
77.72 
33.09 
: 4t.il 

-*.« 
30.95 
0.77*91 0* 
o.iotor 01 


SfiiM t. <01057 
tf 0*19 i. WKHI 
tftiW l. IMCKI 
47*019 C. 

Sftotn o. 

770019 >. »9Um 

coo** jrr* «9 |chi 
W1D.JW.lv l» I 
f*(9«$f»«2. :(• I 
091 X tOKOIfSt 
51690 XX t07( » 
i|W» of loill 
• too *7 lull 

71110 JttfiMUU 

f l*ff 1 

i/vtoie swim ' 


200.00 
-950.00 
255.0 J 
311.00 
-300.29 
935.71 
3*9.71 
-374.29 
935.71 
•070.00 
590.0* 
-997.99 
70.03 

33.14 
30.01 
-2.37 

77.15 
0.30217 02 
0.97*91 00 


229*00 
-929.00 
210.31 
-299.00 
-337.9*: 
290.5* 
-223.** 
-059.4* 
290.9* 
-II* .00 
301.3* 
- 0 * 2 . 2 * 
07.19 
20.13 
20.40 : 
-0.*1 
30.17 
0.252*1 02 
B.M{» 00 


NlOOe««7U01lO<6.O*M0 » 
370*10 1. wm» 
$70019 «• <91*31 
310019 0. <97071 
370019 4. WH« 
170019:0. 1*10*1 
37*019 C. IDKM 
SWttlt. IfffCNR 
97*419 7. («K*t 
'6499* 07. IKfCMI 
99(9.370.1. <• » 
9*Xw.5f*<2» «• I 
49(0 <006*111) 

9(694 XX mil 

5169* 77 . (OSH 

too ox wsti 
tim swns i«it> 

■.1:119511 •:■•'■ ■ 

IWIKilWjai 


2*5.00 

-*24.00 

105.00 

303.01 
-♦79.61 

291.39 
*00.39 

-ii5.*i 

292.39 
-924.00 

4S:S 

41.7* 


■ 1 ; ■: ttflOiOOfOM'l 000.9**901 
314019 0* 190*31 
■ '(99051 ■ 

i«!!£ SKii 

a feS 

9:0* : ;W9ttW ; ' :■ 

9 t. (9CC#» 
19t!09:( ; ; 
374.1. f* » 
(Wife!* < 

w$m « (wootioo , 



*««* n i«||i 

UCntn 

«* Will 
ntno st«m ix*! » 
i i*in 

t/rmn sthw 


'■ «»*93 

r«ffWU*HDie.MMI 
smn «. iifiM 
«*m i. WMl 
<t*«s t. 

5l**lV *. l«fcCH| 

«*nv n iurtrt 
VKmm t. 

$W*«n, l«f(M 

sM-ai* r* i»(<cni 
cm w. 

(• I 
{r* i 

#ffl t ; 

t«$o 

’ fKti n Kill 

fw *« tmn 

«H? fTttii *«*»* 

I : 

tmnn «w$* 


ir*a<^ ir ptiHi 

f : ;y 

*r*«ir* •*« we** 

mktmM mt** 

'/Jfl&w#*: f, 

i«c*i ' 

ttktlH U w«wt 
t#w** **, wt«* 

■ imfti&t+i+.i* * 

try 

*?« I 

M^m w* :, r-y 


'»tt-«r' r ' «*ii» .v : " 
'•mo wm«*n 

fffceS' w»i*» 


i£» : T ; ■ '■ * 


•I m«» 

in»99 

>U!*W 

»«**» 

tn»t* 

1**1? 

tt.»l 

:■• . . *«**/■ 

; (Mi 

W 


~44.W 

40U0C 

443*«0 

*ju.;e 

•»*W 

~ian.cc 

1M.I1 

».» 

D»M 

2S*f* 

. : *).ll : 
11.M 
t.MH « 

).mn m 




Tt**t AfttURil I 

SIUIM. |#C«»I 
ItMIM* l«**| 
ftkllMC. IDCltl 
'?K»W '*» IKKKI 
1 WM V. t*CC*l 
C. INKftl 
MWfl I. t*»KNI 
|M«CHt 
ilMt ttt (Nteto 
MlN.in.1.1 1 » 
Mtn.iw.u t* » 
m» * totems i 

lt(M XX (till 

1*««* W l«f» 
iliu dr iKlti 
v»ite jr«ii will 
i mit 

|/t|H4 »K» 


tM.N 

-m.oo 

-405.00 

205*24 

<m.!t 

l«.W 

•1U.H 

•itUf 

-1MO.OO 

MMI 

14*9* 

tt.il 

40.14 

o.i«>u 
g.MSH 00 


iWUN t, (Mill 
IIMlt t» IMEOSI 
51**10 C. t»M1t 
itUIH I. IWCMfr 

tt«W •• 

5ffc*|0 €. (KfCHI 
ItHlN X* 14fC*l 
Itum V. (MtCMI 
***** XV. 

l 

NtMIM. I** » 
W! X IMIHIU 
J!W» «X (Kill 

tii*t vt mu 
MirtY mu 

VIHO 1TKU ««M» 
I Ixll I 

tmu<> Jtw« 


100.99 

n#»Mu«riKt.rtMi 

iwtix i. i«f*si 
inm », i»i«t 
it»*w c. wm 
iw«w *. »«kw 

tlMW »• IWXI 
lIUMti l*KM» 
57**1’* K. 0*10*1 
liutl) V* IWfCHI 
6*1** ffi t*H»l 


im.m 

1*5.00 

♦TO. 9* 

•Mt.M 

»«» 

•5*0*09 

*o.oo 

214.00 

-11 .12 

•HM1 


•402.05 

».*i 

200*15 

tkll 


W.«l 

•402.05 

•n.w 

•10*0.00 

to»,w 

>24.0* 

♦tl.» 

•m,»i 

»?.*t 

4**27 

t1.» 

10.22 

H.?0 

0.25 

>«.;» 

-5.0* 

M.!t 

20.44 

0.t«M « 

2.000*2 U 

I.IMtf M 

0.220*2 00 

120.01 

120.42 


-07*. 00 

mo 

-*51.00 

v*.t*i 

221.42 

**.1t 

•202.51 

tr.ti 

•510.51 


221**0 

U.U 

-22.51 

*6.11 

•510.51 


•1014.00 




it&V'v'; 


r i 

Ml * tOttOfFSl 
lim n i Kin 
$r«w w mill 

1M « txm . 
omo mi» iuii 
i: iKsti 

imrni jtwjj 


n*n 

**.*t 

»M? 

1*.** 

-T„|4 

e*tW4* oj* 
t.m?f i"j 


■ too ,90 

unn«tTDKiK(.nwi 

twin *. wui 

«t#MN *. IMIOSt 
W«» t. 4KHI 
UuW t, i«rcttt 
i!*«N ». mt«i 
tfMIN C. mrCH» 
titun x* 

SMiiN t. txf cm 

fWU XT* IftCCMt 
NIMt«»U t* ) 
90t*.tTft.2» «• I 
M x lOtftttU 
$IIM XX lull* 

JfOOA W Will 

wit tun 

Ttftt ffUtl t«Stl 

« tost i < 

l*TICL» SWIi < 


105.M 

2.04 

9.Hfl 

144.43 

«.01 

lif.ot 

. IW«81 


■ mm 

fiHW MtUXtlOf C.M** I 
Hunt «. nttii 
ffXMM X. tMfXfl 
wit* c. unti 
«f*M« «• WKN 

fTXAIM 0 * t«CHl 

c* imchi 
itutx x. tutw 

tt*»i» ?. m(ta 
. ♦mu it, mcHi 
Ttl««t 1 » 4 * I* I 
** i 

tut x tottimi 

UVtl JUt fxstt 

IttM IT lull 

»W XT Wtll - 

UHtB twu tun 
imntttuii 


-40.00 

*440*90 

w.ot 

279.09 

J.t* 

*942.91 

W,M 

20.00 

4. At 


t.«l 

-1*41 


*1.44 

70.04 

M.U 

e.fcm* n 

0.222*9 01 

c.IOai* u 


**.oo 

*5.00 

u.«o 

Uft.lt 

: ruu 

*m.« 

ItltW 

192.00 

«.« 

*900.21 


*400.21 

••.TO 

149*29 


-12.21 

5.29 

, -400.21 

• . *99.04 

-040.00 

119*94 

241.54 

-t.n 

*042.04 

19*44 

20.19 

■•‘•V:' 4.99 

9.01 

t*n 

*9.99 

*#•29 

-J.54 

M.tt 

lt.lt 

MM« 91 

«.nw oi 

•*44999*41 







Wm 



mtsi 

«»«N t . UC«D ‘ 
VT**!’* t. IUCCHI 
*4. (Will 
muis t . 1«KM 
itmw »i iwchi 

StMIH V. Utffl 
6MW* rt. l«tV| 


6MM rt. l«tV| 
I * » 

•• u . vtv . 2 . «• * 
Mt i 

SlMtl IH 4I5H 
OK** TV DSD 
Uv XT («Mi 
flftft tTKftS USD 
I UID 

irtifio m « 


I - H.H 
> 4.M 
- 4,00 
wt*oa 

0 . 0t 

1 , » 
03.H0 
OO.W 

1.H 

-•♦.do 

*n.« 

M.w 

I.H 

-'AM 

H.M 

J . I ' W - 1 ! 


1 * 00,00 

♦mi* ♦. i*e*<u 
imw *, w «» 
Murk <» | rt»si 
i»»Mk #, lOiriu 
tiuh >i «#KHl 
* t «* l * t . 

MUIk fc wrjto 
WMk V. DICH* 
M*«» IV. IIICHI 
Hh . Vtf . l . t* i 
Mtk.ttt.i. •* t 
•hi i loi^rtn 
M6U XX USD 

Vtcrt* TV USD 

TM.XV usd •;: 
« m » usd 
I Util 
tmao ir»B» 


00.90 

1.00 

131.00 

0 . 1 * 

*J.|4 

- KM1 
*10. SO 

2 *. 0 * 

Ml 

mi 

- 0**1 

30.0* 

0 . 1101(01 

•.lllrtll 


H.N 

-«««.oo 

-MT.09 
1 * 1.00 
- m.il 
•Mi.ii 
1 * 0.42 
- 10.30 
-Ml.lt 
• 0 * 2.00 
1*2,2* 
-• 45.03 
2V.lt 
•Ml 
-*. 02 : 
• 1.11 
10.03 
4.1HH 01 
J.MW-M 


*3.00 
-002.31 
- 120 . 00 ; 
101*00 
-m.*e 
-♦to.** 
102 . 0 * 

: - 0,*0 

-♦to.** 

•030.0) 
2 * 2 . 0 * 
-•* T.*T 
22 .** 
•!.*3 
-*• J3 
•Ml 
30.0* 
••43021 01 
o.uooe oe 


lifo-i-XA 


It «UH 




CtUfO C901* «M Min 
o * 3t uc . ixfcyna * tut * 


*4*0 tVIEO.MMV MIM 



TEST HO. 5 (COKTIWED): 3.00 INCH TRANSVERSE DISTANCE FROM THE CZXTOtLIHI AT CDfTEl OP PLATE a 
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